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ABSTRACT 

Aims. We study the perturbations of a relatively close third star on a tidally distorted eccentric eclipsing binary. We consider both 
the observational consequences of the variations of the orbital elements and the interactions of the stellar rotation with the orbital 
revolution in the presence of dissipation. We concentrate mainly on the effect of a hypothetical third companion on both the real, and 
the observed apsidal motion period. We investigate how the observed period derived mainly from some variants of the 0-C relates to 
the real apsidal motion period. 

Methods. We carried out both analytical and numerical investigations and give the time variations of the orbital elements of the binary 
both in the dynamical and the observational reference frames. We give the direct analytical form of an eclipsing 0~C affected simul- 
taneously by the mutual tidal forces and the gravitational interactions with a tertiary. We also integrated numerically simultaneously 
the orbital and rotational equations for the possible hierarchical triple stellar system AS Camelopardalis. 

Results. We find that there is a significant domain of the possible hierarchical triple system configurations, where both the dynamical 
and the observational efi'ects tend to measure longer apsidal advance rate than is expected theoretically. This happens when the mutual 
inclination of the close and the wide orbits is large, and the orbital plane of the tertiary almost coincides with the plane of the sky. We 
also obtain new numerical results on the interaction of the orbital evolution and stellar rotation in such triplets. The most important fact 
is that resonances might occur as the stellar rotational rate varies during the dissipation-driven synchronization process, for example 
in the case when the rotational rate of one of the stars reaches the average Keplerian angular velocity of the orbital revolution. 

Key words, methods: analytical - methods: numerical - celestial mechanics - stars: binaries: close - stars: individual: AS Cam 



!>■ 1. Introduction relativistic apsidal motion. There is a small subgroup amongst 

1 : 1 1 1 these eccentric eclipsing binaries which have an additional im- 

p><; In a previous paper dBorkovits et al. | we introduced a stance, as their apsidal advance period is significandy (by 

H . new numerical code which integrated simultaneously the orbital ^^^^ ^^^^ ^^^^^^^ ^ relativistic apsidal motion 

. . ^ , equations of hierarchical triple systems, including the tidal and contribution. It is well-known, that the period of AME in these 

dissipative terms, and the Eulerian eqiiations of stellar rotation. ^ ^^^^^ ^^^^ confirmation or even as challenge 

First we applied the method for the Algol system itself. In this ^j^^ q^^^^^^j ^^^^^^^^ Theory. Unfortunately, these binaries 

triple the inner binary had an almost circular orbit with approx- necessarily have larger separation, so in such systems the ap- 

imately synchronized stellar rotation (another apphcation of an ^^^^j ^^^-^^ -^^ ^^jj^ -^^^ ^^^^^ centuries or even of 

earlier version of the code, not including tidal dissipation was ^^^^^-^ Consequently, in these systems first we have to solve 

g|rnresented for the ternary system IM Aur in| Borkovits et alj ^^^^^^^ ^^curate determination of the apsidal period 

' from a small portion of one revolution of the apsidal line, be- 

In the present study we concentrate on a dynamically less- f^^^ ^an label them as a challenge for the General Relativity 

relaxed scenario, namely when the inner binary has a signifi- xheory 
cant eccentricity, i.e. the system is far from its synchronized and 

circularized state. Perhaps the most important feature of such For our study we chose the eclipsing binary AS 

systems is the apsidal motion eff^ect (AME), i.e. the revolution Camelopardalis, which is a member of an even smaller sub- 

of the orbital axis with a constant period which is determined group of the previously mentioned small group of the eccentric 

mainly by the orbital separation, eccentricity, masses and the in- eclipsing binaries, as this system, together with approximately 

ner mass-distribution of the binai-y members. Nevertheless, sev- six others, shows a significantly lower apsidal motion rate than 

eral other physical processes also force AME. The two most what is calculated from theory. Since this disrepancy was found 

significant ones are the perturbation of a third body, and the for the first time at DI Herculis (Semeniuk 1968), several au- 
thors have investigated this phenome non. A summary of their 
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planations is the perturbations by a third component. The effect 
of the perturbations of a tertiary for the apsidal motion period 
has akeady been investigated only in a few previou s studies 
jKhodvkin & VedenevevI 119971 iKhodvkin et"aLll2004 and ref- 
erences therein). These papers mainly focused on the above- 
mentioned two systems. Furthermore, in our opinion, these ear- 
lier studies have two fundamental disadvantages. First, the third 
body effect and the tidal effect were considered independently, 
and the resultant net apsidal motion period was calculated sim- 
ply as an algebraic sum, which is very far from the reality, as it 
will be shown in the present paper Second, the relation between 
the observed parameters and the physical quantities were not in- 
cluded in the scope of these papers. However, as we show, how 
the observed quantities, which are mainly deduced from some 
variants of the eclipsing O-C curves, relate to the real apsidal 
motion perio d, is a sophisticated problem. Note, as we know, 
IClarell(ll998h was the first to mention this problem, nevertheless, 
in his paper this was not examined in the light of the perturba- 
tions by a third body. 

In this paper we mainly focus on the short term observational 
consequences of the perturbations of a third body for an eccen- 
tric binary. We carry out both analytical and numerical studies. 
We give an analytical form of the time-dependence of the orbital 
elements of the close binary both in the dynamical and the ob- 
servational reference frames up to fifth order in the eccentricity 
and related quantities. We also present the analytical form of the 
O-C diagram of such a binary. We show that the complex varia- 
tions of the orbital elements on a time-scale similar to the tidally 
forced apsidal motion period may result in significant discrep- 
ancies in the shape of the O-C curve from the pure eccentric 
two body case, even without the remarkable real variation of the 
apsidal advance rate. Finally, we caiTy out some longer-time nu- 
merical integration to investigate the variation of the orbital as 
well as the stellar rotational parameters with and without dissi- 
pation. 

It is important to note, that we restrict ourselves only for 
the simultaneous investigation of the third body and the tidally 
forced perturbations in the orbital elements in the frame of 
the classical, Newtonian mechanics. It may seem contradic- 
tory that despite the fact that our purpose is to give some ac- 
ceptable explanations for the anomalously slow apsidal motion 
for those systems, where the relativistic contribution is remark- 
able, we do not consider the relativistic apsidal motion contri- 
bution at all. However, this contradiction can be resolved eas- 
ily, as follows. The previous papers considered the three effects 
(tidal, third body, relativistic) as being independent. If we ac- 
cept this, then our results should only be modified by some 
additive constants, which does not influence our qualitative re- 
sults. Nevertheless, in the following we show that in the sense 
of the tidal and third body terms this is not the case. Similarly, 
we can assume, that neither the third body nor the relativistic 
terms can be considered independently from the tidal contribu- 
tion. Consequently, to correctly consider the varying relativistic 
apsidal motion rate we would have to apply the relativistic for- 
mahsm, which is far beyond the scope of this paper. However, 
in contrast to the tidal term, the relativistic apsidal motion has 
a notably smaller dependence on the eccentricity. This suggests 
that the linear, additive approximation is more realistic in this 
latter case. Consequently, we believe that our calculations give 
significant and well-applicable results. 

We also omit the investigation of the effect of the non- 
synchronized and even non-aligned stellar rotational axes on the 
perturbations of the orbital elements. Theoretically, it can be ex- 
pected that such relatively young early-type eclipsing systems 



such as e.g. DIHero r AS Cam could have non-aligned rotational 
axes (cf.'Zahnll973), which can produce even a reverse d net ap- 
sidal r evol ution. Analytical formulae are given e.g. in IShakural 
(119851) . and lCompanv et al.1 (1 19881). However, fo r the case of AS 
Cam in the thorough discussion iMalonev et al.l (Il989i) showed 
that this solution might be excluded, as the V rot sin ; values de- 
rived from radial velocity measurements of Hilditch ('1972b') 
strongly suggest nearly synchronized rotation. Claret ( 1998) also 
refutes this solution in the case of DI Herculis. We also note that 
although in the presence of a third companion, stellar preces- 
sion could be forced by the misalignment of the orbital planes 
even in the case when perfect synchronization is expected. A 
small amount of amplitude precession of the rotational axes ac- 
tually occured in our numerical integrations presented in Sect. [3] 
Nevertheless, their amplitudes are so small that they could not 
affect the apsidal motion significantly 

In the next section we give the general mathematical form of 
the orbital elements and the O-C curve of an eccentric eclips- 
ing binary when the revolution of the stars are affected by both 
tidal interactions and third-body perturbations. Then in Sect. [3] 
we present several short-time numerical integrations with differ- 
ent initial configurations of the AS Camelopardalis system for 
supporting the analytical results of Sect.|2] and, furthermore, we 
also study the dynamical evolution of the system on a longer 
time-scale, including also dissipative forces. In Sect. |4] we fur- 
ther discuss our results and conclude. Finally, in AppendixlAlwe 
describe our mathematical calculations in details. 



2. Mathematical form of the O-C in a tidally and 
third-body perturbed eccenteric eclipsing binary 

In a previous paper jBorkovits et al.ll2"003h we calculated the ef- 
fect of the third-body perturbations on the moments of the eclips- 
ing minima of such eclipsing binaries which are members of 
close hierarchical triple stellar systems. Here we mainly follow 
the same method described that paper, so we give here only a 
brief summary, except the steps where we substantially modi- 
fied the earlier methods. 

2.1. General considerations and equations of ttie problem 
As is well-known, at the moment of the mid-eclipse 



U X + h 2k7T, 

2 



(1) 



where u is the true longitude measured from the intersection 
of the orbital plane and the plane of sky, and k is an integer 
An exact equality stands only if the binary has a circular or- 
bit, or if the orbit is seen edge-on exactly (for the correct in- 
clination dependence of the occ urrence of the mid-eclipses see 
IGimenez & Garcia-Pelavoll983h . This latter condition is almost 
satisfied in those binaries which are of interest to us now. It is 
known from the textbooks of celestial mechanics, that 



M = — - flcos /, 

PI 

,l/2„-3/2. 



= ti'^^a-'^\l-e^)-^'\l +ecosv)^ -Qcosi, 



(2) 



consequently, the moment of the A^-th primary minimum after 
an epoch can be calculated as 



r-'n n 



Au 



^i/2[l+ecos(M-w)]2j_Pi^^^^.' 

Cl 
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.3/2 



(1 - e2)3/2 



1 P^r^ 
1 H £2 COS I 



1/7 ri / ^ ■ '^"^ 

^ ^ [I + e co5{u - (JL))Y \ C\ 

In the equations above c\ denotes the specific angular momen- 
tum of the inner binary, p\ is the radius vector of the secondary 
with respect to the primary, while the orbital elements have their 
usual meanings. Furthermore, in Eq. (O we applied that the true 
anomaly can be written as v = m - w. Nevertheless, to avoid any 
confusion we emphasize that the angular elements (i.e. u, u, Q, 
are expressed in the "observational" frame of reference, that 
is, its fundamental plane is the plane of the sky, and m, as well as 
<jj is measured from the intersection of the binary's orbital plane 
with that plane, while Q. is measured along the plane of the sky 
from an arbitrary origin. In order to evaluate Eq. (O first we have 
to express the perturbations in the orbital elements with respect 
to u. 

It is well known from the basic works of the three-body prob- 
lem that in the present problem the perturbations in the orbital el- 
ements are effective on three different time-scales. Nevertheless, 
the so-called "short-term", as well as the "long-term" perturba- 
tions can be omitted due to their small amplitude. Strictly speak- 
ing, the second kind of the above perturbations might reach th e 
limit of detectability in some systems (see Bork ovits et al.l2003l) . 
but from our point of view the "apse-node" terms have an ex- 
clusive importance. So, in what follows we concentrate on the 
so-called "apse-node" time-scale perturbative terms. They can 
be divided into two groups according to their different origin in 
Eq. ([3]). First, the "apse-node" time-scale perturbations in the or- 
bital elements e and u) arises also naturally in the formula above 
(as is well-known, there are neither "apse-node"-type, nor secu- 
lar perturbations in the semi-major axis a). We will refer to this 
group in the following as indirect perturbations. Furthermore, 
some other terms which represent low-amplitude, short-period 
perturbations in a, e, o) give large-amplitude "apse-node" terms 
in M due to the multiplication with some of the cosnv terms. 
These are the direct perturbations in the orbital motion. Although 
our calculation of these latter direct perturbations would give 
back the first group too, we found that it is more convenient to 
calculate the two groups in two different ways. 

First, we consider the indirect perturbations. As one can see 
later (e.g. Eqs. 1841 . 1851 ). the variation in both the eccentricity, 
and the argument of periastron during a few revolutions can be 
expressed as 



Aw 



Am, 



(4) 



(5) 



so it is a quite good approximation to carry out the integration 
Eq. ^ first for one revolution treating a, e and w formally as 
constant. In this case taking into account only the first term on 
the right hand side (rhs) we arrive at an analogue of the well- 
known Keplerian equation, which has the following closed solu- 
tion 



Pn^ — 



In 



2 arctan 



2 arctan 



1 - e cos Lo 
1 + e 1 + sin w 



(1 



2)1/2. 



e COS, to 



1 



• cos Cl) 



1 -H e 1 - sin w 



2x1/2 



1 -I- e sin w 
ecos w 



1 - e sin w 



(6) 

for the two types of minima, respectively. (Here P denotes the 
anomalistic or Keplerian period which is considered to be con- 
stant.) Note, that instead of the exact forms above, naturally its 



expansion is used widely (as in this paper), which is as follows, 
up to the fifth order in e: 



P,j, ^ PsE+ — 
2n 



1 . /3 2 1 4. • . 

±-7r + 2ecos oj + —e H — e sin2w 

2 4 8 



/I 1 1 i\ o 5 4 . , 3 , 

+ —e H — e cos 3dL» e sin4w h e cos 5w 

\3 8 / 32 40 

(7) 

where P^ is the sidereal (or eclipsing) period of, for example, the 
first cycle, and E is the cycle-number. Nevertheless, the differ- 
ence of the two quantities (usually denoted by D) is often quoted 
in the literature in its closed form. 

We now formulate the direct perturbations. To do this we 
write e cos v as 



ecosv - e cos w cos M + e sinw sin M, 

- (e cos cl>)o cos u + (e sin a)\) sin u + 

"Ide dco , , 

— cos ix) - e — sin oj \ du cos u + 
du' du' 



(8) 



duj \ , . 
sin w -He — cosdLi dw sinw. 
du' ) 



(9) 



When we concentrate only on the short-period terms in the 
derivatives (i.e. those which are functions of cosnu or sinnM), 
we can take the cos oj and sin w terms out of the integrand, and 
carry out the integrations only for the derivatives, so the "apse- 
node" time-scale direct perturbations in e cos v could be derived 
from 



e cos v'tjii- = cos V 



f {—] dw'-sinv f fe— I du', (10) 



where the subscript „ refers to those terms which contain +u in 
their arguments, and, generally 



e'" cos «Vdir - me'" ' cos nv 



-ne'" ' sin nv 



J»o \ du'l 



du' - 
du'. 



(11) 



Furthermore, the direct perturbations coming from the semi- 
major axis can be calculated as 



~ 2a (1 -Hecosv)2 



f ^du'. 
J du' 



(12) 



The derivatives are as follows 

da 2a^'^ 



dt 



du 



de 
du 
dco 
du 



dQ 
du 



[fie sin v + fi(l+e cos v)] — 



2a^ 



1 



yU (l + e cos v)2 
fl2 (l-e2)2 



[fie sin V -H /t(l + e cos v)] . 



fi (l + e cos 
fl2 (l-e2)2 

fi (1 + e cos 
e cot/ sin M 
1 + e COS V 

fl2 (l-e2)2 



/ COS V + e \ 

fi sin V -H /t cos V + 

\ 1 + ecosv/ 



-fi cos V + /t sin V + 



sin V 



\ + e cos V i 



■/n 



e sin u 



yU ( 1 + e cos v)2 ' sin /(I + e cos v) 



(13) 
,(14) 
- 

(15) 
(16) 
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where /r t,n represent the radial, transversal and normal compo- 
nents of the perturbing force (see later). We applied the follow- 
ing approximation: 



df 
du 



(17) 



Finally, considering the last term on the rhs of Eq. (O, by the 
use of Eqs. ( fTSl ) and ( fT6b it can be seen that all the direct "apse- 
node" terms in e'" cos nv-s which would occur from the /„ force- 
component are cancelled by the opposite-sign term in ^ ^ cos 

and only the indirect term from yU^'^^a^^^^ cos / gives a further 
contribution. 

As a next step we calculate the secular perturbations in the 
above-listed orbital elements. For the calculations we truncated 
the perturbing force at the second order term. In this case the 
force components effective on the close binary are as follows: 



3 Got3 Pi 



8 pI P2 
+(1 -/)^cos(2m' 



[(1 + if cos(2u' - 2u - a) 



• 2m -yS)+2(l - /^)cos(2m - 2Mm)] ,(18) 



3 Gm3 pi 

•^■"2 - 7 2 



(1-/2) cos(2m' - 2u'J + 1^-1 



I Pi 



/r3 = 4 

Pi 

_ 3Gmp\ 
8 pI pi 



(19) 



(20) 



[(1 +lf sin(2M' - 2m - a) 



(1 - /)2sin(2M' -H 2m -yS)-2(l - f) sin(2M - 2Mni)] , (21) 

{2 cos(m - Mm) sin 2(m - Mjjj) sin -I- 
[1 - cos2(m' - M^)] sin(M- Mm)sin2/n,). (22) 



4 pI P2 



We divided the radial component into three parts, relative to 
their different significances. Namely, /i3 is the tidal term, while 
amongst the two three-body terms, fi-2 is formally analogous to 
/r3, i.e. gives a similar secular contribution to the considered per- 
turbations (see later). Furthermore, / stands for the cosine of the 
mutual inclination (/„) of the two orbits, p2 is the radius vector of 
the tertiary, m, m' are the true longitude of the close and the wide 
orbits measured from the plane of the sky, while m^, u'^ are the 
true longitudes of the intersection of the orbital planes measured 
from the plane of the sky along the close and the wide orbital 
planes, respectively (see Fig.[T]). Finally, the contributions of the 
mutual tidal, and the rotational oblateness are as follows: 



T2 = e'^k^^^Rl + '^kf^R' 
\mi 1112 



2 ' 



Grrii 

;,(1.2) 



G;«2 



(23) 
(24) 



where ' are the usual apsidal motion constants, R12 are the 
average radii of the stars, while <y,^ ^ are the rotational angular 
velocities of the star (which are treated as constant here). 

As in the present approximation there are no "apse-node" or 
secular changes in the orbital elements of the tertiary, we consid- 
ered its orbital elements as constant. Although the equations of 
perturbations can be written directly for the above-mentioned or- 
bital elements, it is more convenient to use the equations for the 
orbital elements expressed in the so-called dynamical frame of 




Fig. 1. The spatial configuration of the system. 



reference, in which the fundamental plane is the invariable plane 
of the triple system. The angular orbital elements in the observa- 
tional frame can then be expressed from these by the formulae of 
spherical geometry. To avoid any confusion, the argument of pe- 
riastron, and the longitude of the corresponding ascending node 
in this system are denoted by g and h, respectively (which are 
their usual notations in the perturbation theories). Furthermore, 
the inclination of the orbital plane of the binary with respect to 
the invariable plane is denoted by i[. It can be clearly seen, that 
the relation between the two periastron elements are oj = g + 
(see Fig. Q In the frame of our approximation the studied 
problem is reduced to one degree of freedom. Consequently, first 
we express the variations of all the interesting orbital elements 
in the function of the argument of periastron g of the binary in 
the dynamical reference frame, and then, we give the g{u) func- 
tion. So, the secular parts of the perturbational equations are as 
follows: 



dM 



de 
dh 



1 - ^Qcosi 

ci 

A + B COS 2^ ' 

1 dQ 



■ cos 



dii 



— = tan; I 



A + B cos 2g dg 

eAt sin 2g 
A + B cos2g' 

1 Ani + An2 - Aii2 COS 2g 

cos/i A + Bcos2g 
An2 sin 2g 



A + B cos2g' 
e de 

= rCOtim — , 

1 - d^ 
d<5 Ad + Z?d cos 2g 
dg A + B cos 2g 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 



' Note, that in some studies on apsidal motion for the equivalents of 
Eq.[6]the sign m is used, which is meant as g + h. Naturally, it is correct 
if the perturbative force lies perfectly in the orbital plane of the binary, 
which is fulfilled as far as only the tidal effect of the distorted binary 
members with perpendicular rotation axes are considered, nevertheless, 
in the present situation this is no longer the case. 
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where 6 denotes the direct perturbations in («) while for the 
angular elements in the observational frame of reference we ob- 
tain that 



dQ dh cos Iq - cos ii cos / d/ 1 sin Iq sin h 

dp 1 - cos^ i ' 



dg dg 1 - cos^ / 
d(jj 



1 + 



du„ 



dh dQ 

1 H COS /i COS / 

dg 



1 + 



dh cos /] - cos /() cos / d/ 1 sin /o cos / sin h 



1 - cos^ ; dg I - cos^ / 

A - Ani - A„2 + A, cos2g dQ 



■ cos I, 



dj_ 

dg 



A + B cos 2g dg 
dh sin /o sin ii sin h 
dg sin / 

d/i cos Iq sin ii + cos /i sin Iq cos h 



dg sin / 

In the equations above the tidal contributions are: 



^dr3 = ? 



3 „2 , I „4 



5T2- 



1 + 5^ + „ 
(1 -e2)5 



+ 2';?- 



1 



■(1 -e2)2 



1 + ^^2 + 59^4 + 113^6 

5 16 il 



(1 -e2)5 



1 + 5,2 + 5^4 + 5,6 



(1 - e2)2 

while the third body terms are as follows: 
15 m3 iP" 



Ac 

^di-2 

A, 



8 m 



-(-) 

123 \P'I 



'2N-3/2 



3\ 8 8 64 ' 
Ag(1 

51 31 2 23 4\ 

— e 1 + — e + — e At, 
20 I 51 48 / 



Anl = 



An9 = 



5 



AG(l-e^) 

«2 



2n1/2 



'(1-^2)1/2 

and, finally, 

A = A,3 + Ai.2 + A„i + A 
1 



C2 

El 



2fl5 

+Ag 



3 „2 , 1 „4 



5r2- 



(1 -e2)5 



+ 2?? 



1 



1 



(1 -e2)i/2 



B 



in2 



(1-eO 



2\l/2 . 



/2 - 1(1 



1 



(1 -e2)2 

2y/2 



(31) 



(32) 



(34) 



(35) 

(36) 

(37) 

(38) 
(39) 
(40) 

(41) 

(42) 



+ - 



2 1 + Ci 



5(l-e2)i/2 C2 



C, 



(l-e2)i/2 (l-e2)i/2C2 



Ad = Adr3 + Adr2, 

Bd - Adt, 



(43) 



(44) 

(45) 
(46) 



At the calculation of the formulae above it was also used that 
CCi 



COSii - 



sin;i — 



ccr 

\CxCi\_ 
CCi ' 



(47) 
(48) 



where Ci means the orbital angular momentum of the close bi- 
nary, while C is the same for the whole system. Supposing that 
the rotational angular momenta of the three stars are negligible 
we obtain that 



Ci C2 ^ 
cos li - — + —I, 
C C 

C2 . . 
sinii = --^ sin in,, 



and it is well known, that 
(33) Ci = VGmi2fl(l -e2), 



C2 = 



'«12 

mi2m2 
mm 



ylGm\22,a'i\ — e'2). 



(49) 
(50) 

(51) 
(52) 



2.2. Solution for edge-on binary orbits witli small eccentricity 
variations 

2.2.1 . Closed form solutions for the binary's orbital elements 

For the first time we assume an edge-on binary orbit in the ob- 
servational frame, which is a plausible expectation in the case of 
the relatively wider eccentric eclipsing binaries. Consequently, 
at this stage we omit terms multiplied by cos /. We consider 
Eqs. (|25])-(l32l). One can see that if B > A these equations be- 
come singular at certain directions of the axis. This is exactly th e 
case which defines the so-called Kozai resonance (lKozailll962l) . 
We are interested in such binaries where A > B, i.e. this res- 
onance does not occur. In this case as far as the coefficients at 
the rhs of the equations can be treated as constant, or at least 
their variations are small, all the equations have closed solution, 
which, for B i^Q are as follows: 



1 

U - Uo+ — 



1 



: arctan 



A vrTE2 



1 At / 1 + £'cos2g 
eo - —e — In 



1 -E 
1 +E 
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(55) 
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(57) 
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where 
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where 



(59) 



Eqs. ( 53J and (l55T l are equivalent to the results of I Soderhjelml 
(|1984|) . Nevertheless, in his paper the eccentricity equation was 
not calculated, as the eccentricity was considered as strictly con- 
stant. Furthermore, we stress again, that the last equation for o) is 
calculated in the observational and not in the dynamical frame of 
reference. The first equation reveals that in this case we get the 
following constant angular velocity (in units) for the apsidal 
motion in the dynamical system: 



n = A Vl -£2. 



(60) 



Similarly, the secular terms, i.e. the mean angular velocities of 
the invariant node (/z), and the observable argument of periastron 
iio) are as follows: 



1 



"Ho = 
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while the secular part of the direct perturbations gives 
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In the following we introduce the new variable as 
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It can be seen easily that 
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(63) 
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(73) 

(74) 
(75) 

(76) 
(77) 



(61) respectively. In what follows, we omit the subscript from the 
parameters E, &, because these parameters will always be used 
as constants, with the value calculated at e = eo- 



2.2.2. The analytical form of the apsidal part of the 0-C 

By the use of the Taylorian expansion of Eqs. (l67l i. (fTTT i and ( f72] i 
we obtain for the analytical form of the apsidal part of the 0-C 
as follows: 

y P = y^TT - 2;eo(l + eo) cos[w; + (1 + 04)0] 
+ ^el(\ + eof sin[2w; + (2 + 2'TJ)g] 
+i^-el{l + eo)' cos[3wS + (3 + 31/)g] 



1 



1 
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1 
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+ — eg£sin2g + (9[(e,£)4]. 



(78) 



From now on we suppose, that (9(e) - 0(E) = <9(fi), or more 
^ ^ generally, (9(e) = (9 ( ^ ) . Furthermore, 
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Anl + An2 A 
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(79) 
(80) 
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the latter is the period ratio of h cos ii (or in the present approxi- 
mation. Mm) and g. This term gives the relative difference of the 
speed of the apsidal advance in the observational and the dynam- 
ical frames of reference. Furthermore, note, that in this approxi- 
mation the direct terms give only a small contribution to the O-C 
(the last term of Eq. 1781 ). Nevertheless, the secular part of such 
perturbations are added to the observed eclipsing period. Finally, 
j + I for the two different types of minima. 

We consider the two extreme cases. First, if the two orbits 
are coplanar (i.e. = 1), then the angular velocity of the apsidal 
motion is 



(1 + ton = A,2 + Ar: 



(81) 



Furthermore, as fi = 0, (i.e. Ae = 0, which is true as far as we do 
not consider the octuple term in the perturbing force), Eq. jTSl l 
reduces for its usual form, apart from the different period given 
by Eq. dSTT l. Second, in the case of two perpendicular orbits (i.e. 
I - Q), U diminishes, as well as fi = £ occurs, and, conse- 
quently, Eq. ( iTSl l also becomes somewhat simpler Moreover, the 
more important feature is that these are generally the only two 
cases, when the O-C curve has only one fundamental period. 
Finally, in these two extreme cases 



(82) 



so the results above are rigorously correct for not only edge-on 
visible orbits, as far as the orbital eccentricity can be considered 
as constant on the rhs of the perturbation equations (l25T l-(l33Tl. 

2.3. Solution for the general case 

2.3.1 . Formulae for e, g and h - results and discussion on 
the dynamical apsidal motion period and on the nodal 
regression 

One can see that the assumption of the constant eccentricity re- 
mains plausible only if B << A, or if At ~ 0, which may happen 
in two different ways. The trivial case, when no third body exists 
in the system, or at least its influence can be disregarded with 
respect to the tidal forces, or the other possibility is, that the 
two orbits are nearly coplanar. (Note, that this latter case usu- 
ally also could satisfy the B « A condition, as in this case B 
is in the order of e^Ac-) Nevertheless, in the really interesting 
systems neither of these conditions are fulfilled, so we have to 
solve the equations above in several iteration steps. In order to 
do this we used the Taylorian expansions of Eqs. ( |25l l, ( |26b with 
respect to e. Our calculation is listed in Appendix [A] Here we 
give only the final forms up to the third order in the inner eccen- 
tricity, together with such assumption that the other quantities 
given above are also in the first order of eccentricity. So, in this 
approximation the modified angular velocity of the apsidal ad- 
vance, the inner eccentricity, and the argument of periastron in 
the dynamical frame (g) become 



n 



1 , 1 
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(83) 
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(87) 
(88) 
(89) 



Before continuing with the effects of those quantities which 
relate to the observational system, we discuss our results for the 
dynamical apsidal advance rate. First, we consider the instanta- 
neous angular velocity of the apsidal line, i.e IT = VA^ - B^. As 
we are concentrating on hierarchical triple systems, where the to- 
tal angular momentum is highly concentrated in the wider orbit, 
we omit the terms in A and B which are multiplied by Ci /C2. In 
this case we can easily have that up to = 63°43 (i.e. > 0.2) 
A > A,3. So, when 



A,-2 + Anl + An2 > Aj - An2 



(90) 



then the third body affected instantaneous apsidal angular ve- 
locity is surely larger, and, consequently, the apsidal motion pe- 
riod is shorter than in the only tidally perturbed case. Note, that 
Eq. (|90] l is also the zero order condition for the occuiTence of 
the so-called Kozai resonance in a mass-point three-body model. 
This gives < 39?23 for e = 0. We now concentrate on the aver- 
age angular velocity, 11*. There is only one quantity in the 11/11* 
ratio, Eq. ( [83] l which can be negative, namely eo- Consequently, 
all of the other terms would produce a shorter average period 
than the instantaneous one. According to Eq. (|87] |. there is a first 
order term in eo, \S> cos 2gQ, which may give the largest contribu- 
tion in the whole 11/11* ratio. One can see that we could expect 
a longer average apsidal motion period than the instantaneous 
one, when cos 2gQ =s - 1 . This happens when the argument of 
periastron is ^0 ~ +90° in the dynamical system in the moment 
of the calculation of the orbital element (i.e. at the time of obser- 
vation). Note that as one can see from e.g. Eq ( |54l i, for g^ ± 90° 
the eccentricity takes its maximum value, and as is well-known, 
the larger the eccentricity the faster the tidally-forced apsidal ad- 
vance speed, this result is not an unexpected one. 
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For the sake of completeness, we give our result for the dy- 
namical nodal regression, although up to third order in e, E it is 
identical with the Taylorian of Eq. 



2 COS i\ n* 



where 



1 1 



Hq = -f^^- Mm + An2 1 1 + \e\ + \eQ^& 
n*cosii I \ 2 / 4 de 



(91) 



(92) 



We emphasize again, that this term gives further, negative con- 
tribution to the observable apsidal motion period, through the 
/icos/i expression, which is independent of the observable in- 
clination of the system. 

At this point we refer to the paper of iHegediis & Nuspll 
(Il986l) . In their work they investigate the possible observational 
effects of nodal motion forced by inclined stellar rotational axes. 
Although their treatment is entirely coiTect, they unfortunately 
denoted the dynamical node (which is h in the present paper) 
by Q, which usually means the longitude of the node in the sky 
(as used in the present work) in the observational astrophysics. 
From this notati on some misinterpretations occured later So, 
when iKhodvkinI 0989) reacts to the previously mentioned pa- 
per, and states that the orbital plane precession of an eclips- 
ing binary is unable to si gnificantly distort the observed apsi- 
dal motion rate (see also iKhahuUin et anil991h . he is right in 
the sense of the observed node (see the dQ cos / term in d<y), 
but Hegediis & Nuspl (1986) consider the dynamical node (the 
contribution of dh cos i\ in dw). We note also, that in the inclined 
rotation case i\ is small, whi ch justifies the omission o f the cos i\ 
multipUcator in the work of lHegediis & NuspH ( Il986l) . 

2.3.2. Non edge-on orbits: further observational effects from 
nodal motion (w, u^) 

We know continue with allowing non-exactly edge-on orbits 
(i.e. cos / + 0). We then have to take into account the dQcos/ 
terms, too. According to Eq. d3TT i. as far as we omit the small am- 
plitude term di'i this can be written as a function of d/z cos i\ . The 
detailed calculations together with the results up to fifth order in 
e and E, S are given in AppendixlAlCsee Eqs. ||A.114| - ||A.120| ). 
Here we list the form of the result up to the fourth order in the 
secular term, and to the third order in the trigonometric ones. 
Namely, 



to^ 2 + U„ 



where g is given by ( l85l l. while 

«m = (Mm)S + Uog - i^(l + Co) sin 2^ + Ci sinih; + HoQ) 
+ ^C2 sm(2h; + IHqQ) + ics sinO/jJ + 2>Ho§), (94) 

where 

£/o = ~ |a„i + A„2 (l + ^/ij + ^^o^fij (1 + Co). (95) 

The Co-3 terms are coming from the Taylorian expansion of 
(1 - cos^ ' which arises in the equation for dflcos/. These 
are trigonometric functions of /o and ii. They are listed only in 
AppendixlAl Eqs. dXseb - dATTjl ), ( IA.121I ). ( IA.122b . In our cal- 
culations we consider C„ as it would be «-th order in e. 



Here we would like to stress that in this case a further im- 
portant difference to those mentioned earlier is that there will 
no longer be one true period for the apsidal advance in the ob- 
servational system, now we can speak only about only quasi- 
periodicity. (This latter statement applies even better for higher 
orders where the different linear combinations of the angular ve- 
locities of g and h appear) One can see this clearly in Fig. |2] 
where we plotted the variations of the eccentricity and the ob- 
servable argument of periastron (w) in a hypothetical eccen- 
tric triple system for two different initial mutual inclinations 
(i'm = 20°, 60°) of the close and the wide orbital planes. (The 
initial parameters are listed in Tables [T] and |2]) In these figures 
we connected the lo - 2n values with the corresponding eccen- 
tricity values by dashed lines. Here it is clearly visible, that the 
periods are not equal to the corresponding e and o). It can be 
seen also that in the case of /„, - 20° (left panel) the eccentricity 
obtained from numerical integration (see next section) signifi- 
cantly departs from the analytical value, while the argument of 
periastron (w) shows a good fit, at least as far as one omits the 
small cumulative error in the period. On the contrary in the case 
of = 60° (right panel) eccentricity curves produce very sat- 
isfactory fits, while the analytical u) departs from the numerical 
curve suddenly after approx. 1000 years. In the first case, the ex- 
planation can be found in the small mutual inclination, in which 
case, as was mentioned earlier, the higher order terms in the per- 
turbative forces are also important, while in the second case the 
departure arises from the fact, that due to the large amplitude 
orbital precession the cos term becomes so large that our ap- 
proximation will no longer be satisfactory. 

2.3.3. Direct perturbations in the orbital motion 

We also calculate the direct terms for higher accuracy. This is 
slightly problematic, as due to the very strong e-dependence, the 
derivatives of Ad could increase to very large values already for 
medium eccentricities. Fortunately, if the time differences of the 
two types of minima are used instead of the usual O-C function, 
these direct terms will fall out, as at the same time they have 
the same value for primary and secondary minima. Nevertheless, 
we use them in the following, considering the quantities derived 
from the derivatives of Aj as first order, and those from as 
second order (in our present sample configuration of AS Cam 
this can be done approx. up to e' - 0.3 for perpendicular orbits). 
The results up to fifth order can be found in Eqs. ( |A.70l )-( |A.75t , 
while up to third order are listed below: 
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Fig. 2. The variation of the eccentricity and the observable argument of periastron in the case of the AS Camelopardalis due to the 
perturbations of a third body. The initial mutual inclinations are = 20° (left panel), and = 60° (right panel). The dashed lines 
connect the o) - In values with the corresponding eccentricity, illustrating that the period of the e-cycles and the observable apsidal 
motion is different. (See text for an explanation of the discrepancy between the numerical and analytical results.) 
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2.3.4. Generalized form of the 0-C 

Finally, now we report the generalized form of the O-C curve, 
which is as follows: 

y O - C ^ y,00 COS[w*, + (1 + Uq)0-\ 

+y2oosin[2w,* +(2 + 2t/o)^] 
+ y30ocos[3wo + (3 + 3f/())^] 
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+V1-2-1 cos[w; - hi - (1 - Uq + Ho)g] 
+V1-21 cos[w,* + /i; - (1 - Uo - Ho)g] 



+V12-1 cos[w,'; - 111 + (3 + Uo - Ho)g] 

+Vm cos[w*) + hi + (3 + Uo + Hq)^] 
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+V02Q sin 2g 
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+Vo(n sm[hl + HqQ] 
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Furthermore, 



(121) 
(122) 

(123) 



A more detailed result up to fifth order containing 102 trigono- 
metric terms is also listed in AppendixlAl Eqs. ( IA.123I )-( IA.225I ). 
Here we note, that the indices refer to the multiplicators of <y, g, 
and h in the trigonometric terms, respectively. We separated the 
different terms by blank lines. All four groups besides the pertur- 
bations in e and g, give the additional contribution of the preces- 
sion of the orbital plane with respect to the observer^ The first, 
and second groups give the nodal contribution to the observable 
argument of periastron (w), while the third one comes both di- 
rectly from the Q cos /, and the direct orbiatl motion terms in ii. 
Finally, the fourth one also contains the contribution of Q cos /. 
These latter two terms, naturally have e-independent parts, con- 
sequently they do not disappear even in the case of a circular in- 
ner orbit. Nevertheless, as was mentioned above, due to the cos / 
multiplicator, their significance is very limited in the case of the 
relatively wider eccentric eclipsing binaries. Similarly, due to 
this latter condition, the contribution of the second group is also 
very minimal, as these terms arise exactly from the same reason 
as the terms of the fourth groups, manifesting the same preces- 
sion indirectly, through the variation of o). On the contrary, the 
first group contains such terms of the nodal motion which re- 
main effective even in the edge-on cos / - case, too. These are 
coming from the h cos /] term as was discussed in Sects. 12.2.21 
and lzfn 

To illustrate the net effect to the O-C we plotted in Fig. [3] 
the O-C diagram of AS Camelopardalis for the same two con- 
figurations for which the variations of e and oj were shown in 
Fig. |2] The upper panels represent the numerical results derived 
directly from the numerical integrations as well as analytical re- 
sults calculated according to the higher order formula listed in 
Appendix |A] Eq. (|A.123t . The lower panels represent the nu- 
merical curve and the "unperturbed" (i.e. only two-body tidal 
distortion is present, as usual) theoretical curve where both the 



^ We do not use intentionally the usual phrase 'nodal regression', as 
we concentrate on the observational frame of reference, where we can 
observe even nodal progression (see Sect.[3]l. 



constant eccentricity and apsidal motion period were deduced 
from the "measured" eo - e(fo) value. We note, that the signifi- 
cant departure of the numerical and the analytical curves which 
occurs in the - 60° case after approximately 700 years (upper 
right panel) should arise from the fact that due to the variation 
of the observable inclination (;) for this time our fundamental 
principle, i.e. Eq. ([TJ ceases to be correct. Nevertheless, at this 
value of /, our system is already no longer an eclipsing one, con- 
sequently, this situation is out of our scope. 

2.4. Discussion on observational effects 

When only a small fraction of the O-C curve is observed, several 
serious problems occur in the sense of the determination of the 
apsidal motion period with respectable accuracy. To demonstrate 
this we consider only the Vioo and Vi-20 terms (we omit Vno be- 
cause we concentrate mainly on such configurations where the 
mutual inclination of the inner and the outer orbit is large, and 
consequently, we suppose that S ~ E). Furthermore, when the 
period of the nodal motion is significantly longer than the tidally 
induced apsidal motion (as in the case of AS Cam, but not neces- 
sarily in e.g. DI Her), then well as h can be considered as 
constant. Then instead of the usual O-C we apply the difference 
function of the primary and secondary minima, as 



-eo [-2 cos[wo +go + n*(M - uq)] 
n * 

+E cos[wq -So- n*(M - Mo)]} . 



(124) 



As far as the nodal motion is neglected it can easily be seen that 



furthermore, according to its definition 

00 = ^0 - ^£'sin2go + 
so. 



w,') + Qq ~ (A)Q- -E sin 2^0, 
1 

(^0 ~ Go ~ "^0 ~ 2^0 + sin 2go. 



(125) 

(126) 

(127) 
(128) 



We concentrate continuously on perpendicular orbital planes 
configuration. In this case, taking into account that the binary 
is an eclipsing one, i.e. its observable inclination is close to 90°, 
the plane of the outer orbit should lie either close to the plane 
of the sky, or perpendicular to both the planes of the inner orbit 
and the sky. In the previous case (which is the more interesting 
one, because it makes it possible not to observe light-time effect) 
it is easy to see that w s; g + [0, 1] x tt, while in the latter one 
OJ ss g ± 7r/2 (we also used the fact that in these hierarchical sys- 
tems the outer plane is close to the invariable plane). In the two 
different situations we get that 

Ao;r ~ -2— eoflH — £ cos 2(^0 £|cos [wo + n*(M - mo)] , 

;r \ 2 2 / 

(129) 

A±;r/2 = -2^eo^l - ^£cos2wo + ^sjcos [wo + n*(M - Mo)] 

(130) 

The subscript of A refers to the approximate value of m^ and 
we omitted the sin2^o terms inside the arguments. Eq. (11291 ) 



T. Borkovits et al.: Tidal and rotational effects in the perturbations of hierarchical triple stellar systems 



11 






' ' MIN| (num.) 




MINii (num.) 




MIN| (anal.) - : 




MIN|| (anal.) 







1 500 2000 
Years 




1 500 2000 
Years 





MIN| (num.)' 




MINii (num.) 


- 


^m. MINi (anal.) - : 




# % MINjiianal.) 











1 500 2000 
Years 




1 000 1 500 2000 

Years 



Fig. 3. The numerically (AS2 and AS3 runs - see next section) and analytically generated O-C curve in the case of the AS 
Camelopardalis with the perturbations of a third body. The initial mutual inclinations are - 20° (left panel), and /„, = 60° 
(right panel). Upper panels show the numerically generated O-Cs, as well as the analytical ones calculated according to the for- 
mulae of the present paper. Lower panels demonstrate the difference between the numerical (i.e. "observed") curves, and the only 
tidally forced apsidal motion produced O-C. 
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Fig. 4. As in Fig. [3] but the initial mutual inclinations are /,„ x 90° for both panels (AS4 and AS4b runs - see next section). The 
only difference is that in the left panel the orbital plane of the tertiary is almost perpendicular to the plane of the sky, while in the 
right panel this nearly coincides with it. Note that at this second configuration the amplitude of the O-C is definitely smaller than in 
the first case. 



shows a very important result. In the situation when a third body 
revolves around the eclipsing binary in a nearly perpendicular 
plane which lies close to the plane of the sky, the amplitude of 
the A curve - which is used several times for the calculation 
of the period of the apsidal motion - could be lower than it is 
expected (by the multiplicator > 6 > 1 - £), consequently, 
the numerical fitting will result in a smaller angular velocity, 
i.e. longer period. This is demonstrated in Figs. |4] and |5] where 
two numerically integrated O-C curves are shown (together with 
the corresponding analytically calculated ones). There is no dif- 



ference in the physical configuration of the entire triple system, 
only its orientation rotated with 90° along the inner orbital plane, 
i.e. the outer orbital plane from an almost perpendicular position 
became nearly coincident with the plane of the sky. Moreover, 
besides the smaller amplitude. Fig. |5] reveals a more important 
result. As one can see, when the orbital plane of the tertiary is 
close to the plane of the sky, large almost horizontal regions can 
be found in the difference curve. (This means that in these inter- 
vals the primary and the secondary minima vary almost in the 
same manner.) The mathematical cause can be understood from 
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Fig. 5. The difference of the O-C of the primary and the sec- 
ondary minima for the O-Cs plotted in Fig. H] We shifted the 
i' = 178° curve with approx. 249 years in order to coincide with 
the "unperturbed" curve of the two cases. 



the third order approximation. In this case the Vioo, Vi±20 and 
the V'3()o terms should be counted. Limiting ourselves only for 
the third-body-in-the-sky case, then 



Ao,;r 



-2—eo 

7T 



1 - -£+ -£cos2wo + 
2 2 8 



1 , 1 , 1 

--£ cos2<x)o + 770 cos4wo + o^i^' 
4 16 8 



cos^ 



1 



1 



-Ai£+-e^|cos3^ 



(131) 



The flat extrema occur when Q - Q + kn. Around these inter- 
vals the two cosine terms have opposite signs. This also happens 
in the "unperturbed" case of course. However, a significant dif- 
ference is that, whilst we considered both A \ and & as formally 
having the same order than e, both could be somewhat larger 
than e, and in some cases they can reach almost the order of 
e'^^. Consequently, the superposition of the opposite sign cos 3^ 
curve onto the cos Q may cause further flattening almost in the 
order of e. These areas might cover even the 40-50% of the total 
curve. It is trivial, that if this system were observed within this 
time interval, the apsidal motion period would be found signifi- 
cantly longer than the theoretically expected value. On the other 
side, there is no region, where the slope of the difference curve is 
significantly larger than the "unperturbed" reference curve. The 
consequences of the fact above will be discussed in Sect. |4] 

We now study numerically the eclipsing system 
AS Camelopardalis, where only about a 10% or less of 
the total period is covered by observations. 



Table 1. Fixed initial parameters of the binary as well as of the 
tertiary for every run. The k^'^ co nstants are taken from the ta- 
bles of IClaret & Gimene3 (1199 lb . The other parameters for the 
binary ( with the exception of Q which has an arbitrary value) 
are from lKhodvkin & Vedenevevl(ll997h . The parameters of the 
third companion (prim ed quantities) are partly from our earlier 
light-time solution (see lBorkovitsll200"3h . 



a 
e 
■f 



n 



17.195Ro 
0.17 

50000HJD 

50000.82 

50001.7 

88°78 
130° 



a' 736.98Rq 

e' 0.41 

r' 52 085HJD 



mi 



3.3Mo 
l.lMo 
2.6OR0 
0.0049 



kf 0.0011 
^1 0.0 
0.0001 



m2 2.5Mo 

R2 I.96R0 
kf 0.0038 
kf 0.0008 
/I2 0.0 
0.0001 



a: The three values refer to departure from periastron, instantaneous- 
angular-velocity-equal-to-average, and apastron positions, respectively. 



both in the observational (w, u' , Q', /'), and the dynamical {g, 
h, i'l, h' , 12) coordinate systems. Note, our input parameters for 
the integrator are only the observational angular elements, and 
the others are only derived quantities. All columns of this latter 
table refer to four individual runs. One of them is the dissipa- 
tionless case, starting the binary from its periastron, and three 
of them with dissipation (Ai = A2 = 10"'*), with initial posi- 
tion of the binary in its periastron, apastron, and at true anomaly 
V ^ 104°, where the instantaneous orbital angular velocity in the 
present configurations is nearly equal to the averaged one (as we 
used the option of our code which makes the angular velocity 
vectors of the stellar rotation for both stars equally long and par- 
allel to the initial instantaneous orbital angular velocity vector 
in an iterative way, it means that the initial rotation of the stars 
were synchronized in these three diff'erent ways.) We emphasize 
again, that as we were not interested in it exclusively dynamical 
evolution, but for observational consequences too, we carried out 
some runs where the dynamics of the triple were the same, but 
its orientation with respect to the Earth were different (e.g. AS3, 
AS4 and AS3b, AS4b runs). 



3. Numerical studies 

We carried out sever al sets of integra t ions w ith the numeric in- 
tegrator described in iBorkovits et al.l (12004 ). These runs partly 
serve as numerical support for the analytical calculations de- 
scribed in the previous section, and partly serve as study on 
the dynamical evolution of eccentric hierarchical triple systems 
with and without dissipation. The in itial p arameters of the binary 
were taken from lKhodvkin & Vedienevevl(ll997t) . with the excep- 
tion of the i nner structure constants wh ich were set according to 
the tables of lClaret & Gimened(ll99lb . The physical parameters 
of the three stars as well as those initial orbital elements which 
are the same in the different runs are listed in TablelT] while those 
orbital elements which differ in the individual runs are listed in 
Table |2] In this latter table we give the angular orbital elements 



3.1. Non-dissipative runs 

The variation of the orbital elements during the integrations can 
be seen in Figs.|6l|9] These figures represent the dissipationless, 
departed from periastron runs. The left panels cover a time in- 
terval of 100 years, which is naturally briefer than a moment 
from dynamical point of view, but this is the time interval (in 
the best scenarios) which can usually be reached for the present 
studies. The numerically generated O-C curve, as well as the 
variations of the orbital elements e and a> of the AS2-AS4b runs 
were also used in the previous section as an illustration for the- 
oretical thinkings (see Figs. |2]-[5]l. Additional integrations were 
performed for the close binary without the tertiary component 
both in the quasi-synchronized rotator, and the precessing sec- 
ondary component case (this latter is not presented here). 
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Fig. 6. The variation of both the dynamical and the observational 
orbital elements of the binary during 100 and approx. 2800 years 
(1 million days). Semi-major axis is given in Rq, while the angu- 
lar elements in degrees. The mutual inclination is i'm ~ 0?8. See 
Col. ASl of Table|2]for the initial values. 




Fig. 7. The variation of both the dynamical and the observational 
orbital elements of the binary during 100 and approx. 2800 years 
(1 million days). Semi-major axis is given in Rq, while the angu- 
lar elements in degrees. The mutual inclination is » 20°. See 
Col. AS2 of Table|2]for the initial values. 




Fig. 8. The variation of both the dynamical and the observational 
orbital elements of the binary during 100 and approx. 2800 years 
(1 million days). Semi-major axis is given in Rq, while the angu- 
lar elements in degrees. The mutual inclination is » 60°. See 
Col. AS3 of Table|2]for the initial values. 



Table 2. Initial orbital elements, and some derived quantities, 
which were different at individual integration runs. 





ASl 


AS 2 


AS 3 


AS 3b 


AS4 


AS4b 


8 


45=0 
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145? 1 
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352°93 


350:75 


351:00 


172?22 


351? 00 


a>' 


261° 


261° 


261° 


171° 


261° 


351° 


h' 


179?87 


267° 95 


269:89 


1 80:00 


89? 27 


0?00 


n' 


130° 


110° 


70° 


130° 


220° 


130° 




o;i5 


3°75 


10:27 


10:26 


13? 16 


13? 16 


i' 


88°0 


88':o 


88:0 


28?78 


88?0 


178?7 


im 


a78 


20:01 


59:98 


60:00 


89?96 


89? 92 




Fig. 9. The variation of both the dynamical and the observational 
orbital elements of the binary during 100 and approx. 2800 years 
(1 million days). Semi-major axis is given in Rq, while the angu- 
lar elements in degrees. The mutual inclination is » 90°. See 
Col. AS4 of Table|2]for the initial values. 



In the quasi-synchronized binary case U - 381 800 d » 
1 046 y was found for the period of the apsidal motion without 
the relativistic contribution. This yields Wci = ^ci - 34?4/100y 
for the classical part of the apsidal advanc e rate, which is i n good 
agreement with the theoretical value of iMaloiiev et al.1 (Il989h 

= 35°8 ±5°8/100y. 

In Fig. [To] we plotted the variation of the apsidal line both in 
the observational (w), and in the dynamical (g) frame of refer- 
ences for five different configurations of the system. Comparing 
the four third-body perturbed runs with the quasi-synchronized 
binary configuration, one can see that there are no cases where 
the apsidal motion period would be significantly longer than in 
the reference case. However, on the contrary, for the low mutual 
inclination cases (ASl, AS2), the apsidal motion period is found 
to be remarkably smaller. This is in good coiTespondance with 
our discussion on (11*)"' in Sect. 12.3.11 Consequently, we can 
conclude again, that the observational effects due to the specific 
spatial configuration of the orbital plane of the perturbing third 
star should play a more important role in the explanation of the 
anomalously slow apsidal motion than the dynamical perturba- 
tions themselves. 
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Fig. 10. The variation of the argument of periastron in the observational (left panel) and the dynamical (right panel) frame of 
reference. In the right panel we shifted the individual curves in time to start the ^-curves from the same initial value. 
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Fig. 11. The variation of the orbital elements of the binary in the dissipative case. The initial parameters are listed in Table |2] 
Col. AS4. Left: The stellar rotation was initially synchronized in the periastron. Medium: The stellar rotation was synchronized to 
the average value of the Keplerian angular velocity. Right: The stellar rotation was synchronized in the apastron. 



3.2. Dissipative runs 

Considering the dissipative runs, in Fig.[TT]we give the variation 
of the orbital elements in the case of the AS4 configuration for all 
the three initial synchronizations, i.e. when the synchronization 
was done to the periastron (left panel), the average value of the 
orbital angular velocity (medium panel), and the apastron value 
(right panel). The following differences can be observed. 

The shrinkage of the binary orbit is more than three 
times faster during the total integration time in the apastron- 
synchronized case. For a qualitative explanation we recall here 
the simplified dissipative pa rt of the eq uation of motion of the bi- 
nary from Eq. (23) of Borkovits et al.l (,2004.) . According to this 



dis 



3 — 2~^i 

Pi 



Pi X Pi 



(132) 



where pi is the Jacobian vector directed from the primary to the 
secondary, and /?, is the difference vector of the rotational an- 
gular velocity vector of the i-th star (/ = 1 , 2) and the orbital 
angular velocity vector As one can see easily, in the apastron- 
synchronized case the constant ■ p, x pi terms give a continu- 



ously negative transversal force component (which has a maxi- 
mum amplitude at periastron, and almost zero only around the 
apastron). 

The amplitude of the e-cycles is also somewhat larger in 
these latter cases, and the period of the apsidal advance is longer. 
Nevertheless, the difference in the period between the periastron- 
synchronized and apastron-synchronized case is about 15%, so 
it remains clearly under the amount of the observed discrepancy. 
Furthermore, we note that, although the orbital shrinking as well 
as the synchronization of the angular rotation are evident on the 
time-scale of the present integration, it can be clearly seen that 
there is no evidence for the circularization of the inner orbit. This 
suggests that the presence of a not so distant third body may ba- 
sically modify the dissipative circularization process. Naturally, 
further investigations are needed in this question. 

In contradiction to these smaller differences in the behaviour 
of the above-mentioned orbital elements, the angular elements 
(both in the observational and the dynamical frames) show com- 
pletely different variations in these latter two cases than in the 
periastron-synchronized one. The reason can be found in the spe- 
cial initial configuration, i.e. in the (almost) perpendicular posi- 
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tion of the inner and outer orbital plane. In this position the dif- 
ference in the angular momentum stored initially in the stellar 
rotation was able to change the direction of the orbital preces- 
sion. 

Finally, we consider the stellar rotation affected by the dis- 
sipation. As is well known, the phase space of the rotation of 
oblate bodies usually contains large chaotic regions (in the sense 
of our Solar system see e.g. Laskar & Robutel 1993). In Fig. [12] 
we show some interesting resonances. We found in several dif- 
ferent runs, that when the angular velocity of the rotation of one 
of the stars becomes temporarily equal to the average orbital an- 
gular velocity of the binary, typical resonance phenomena occur, 
i.e. the amplitude of the stellar precession suddenly increases, or 
even some fluctuations arise in the stellar rotational rate, and, 
consequently, this can manifest even in some similar fluctua- 
tions in the binary's semi-major axis, as happens in our AS3b 
run. Detailed investigations of such resonance phenomena may 
be the subject of future studies. 

4. Discussion and conclusion 

In this paper we gave analytically the net effect of the perturba- 
tions of a third companion and the mutual tidal and rotational 
distortion of the binary members for the orbital elements of an 
eccentric binary star. We investigated primarily how the presence 
of the third body can affect the O-C curve, which is the most 
important source of the apsidal motion information in relativis- 
tic eccentric eclipsing binaries. For this purpose we have cho- 
sen one particular eclipsing system, namely AS Camelopardalis. 
Nevertheless, naturally, most of our results are valid in general, 
independently of this specific system, although we stress, that in 
the case of other systems, or even in the same system with dif- 
ferent third body parameters, the magnitude of some parameters 
can differ from our assumptions. 

Our work on this topic is not the first, but one of its main 
novelties is that we calculated simultaneously the effect of the 
tidal distortion and the tertiary perturbations. There is a strong 
interaction between the two physical processes. Here we refer 
to the strong dependence of the tidal forces on the eccentricity. 
Consequently, the third body forced eccentricity variation influ- 
ences notably the tidal phenomena, and besides the pure, third 
star forced apsidal motion rate variation, this effect produces fur- 
ther remarkable change in the apsidal motion speed. At this point 
we emphasize again the necessity of the common treatment of 
the two different physical processes, as their simultaneous pres- 
ence can basically change the evolution of a close binary sys- 
tem. We mainly refer to the so-called Kozai res onance, but fur- 
ther co nsideration of this topic can be found in lEggleton et alj 
(1199 8'). Note that from our Eq. ( |26] l naturally almost the same 
condition for the ti dal prevention of the Ko zai resonance arisefl 
which was given in lKhodvkin et alj (|2004|) by Eq. (38). 

The other novelty, as we already emphasized, is the calcu- 
lation the mathematical shape of the O-C curve, and we thor- 
oughly investigated the reliability of the data which can be ob- 
tained from such a distorted O-C when only a small portion of 
the total apsidal motion period is covered. 



' We said "almost", because it is not clear for us why the authors use 
in their analytical stability investigation such orbital elements which 
refer to the observable system. It should be clear that there is no con- 
nection between the physical instability in the system and its orientation 
with respect to the Earth. So in their Eq. (37) which defines the stability 
criterium they should use (according to their notation) (f> (this is -g in 
the present paper), instead of o). 



We found that there is a significant domain of the possible hi- 
erarchical triple body system configurations where both the dy- 
namical and the observational effects appear corresponding the 
tendency to measure a longer, and even remarkably longer apsi- 
dal advance rate, than is theoretically expected according to the 
known physical parameters, and the measured eccentricity and 
visible argument of periastron. This happens when the mutual 
inclination of the close and the wide orbits is large, i.e. the orbits 
are nearly perpendicular to each other, and, furthermore, the or- 
bital plane of the tertiary almost coincides with the plane of the 
sky. As the observable inclination of the anomalously slow apsi- 
dal motion eclipsing binaries is necessarily close to / = 90°, this 
means that almost every second near a perpendicular triple sys- 
tem might belong to the ideal case. Note, that the first condition 
which mainly (but not exclusively) comes from the dynamical 
calculations is in accordance with the results of the previously 
mentioned studies. However, we have shown that it has a notably 
lower significance alone as stated previously. So, we concluded, 
that the observational effects can play a remarkably more im- 
portant role in the detection of anomalously slow apsidal motion 
rate than the dynamical ones. 

This result is all the more impressive, because it offers a nat- 
ural answer for a further serious question: why have we not ob- 
served these third companions yet? We consider this question in 
detail. Several direct and indirect methods can be used now to 
detect a not so distant tertiary in a binary system. A relatively 
detailed, but not complete listing of these methods is found in 
iPribulla & Rucinskil ( |2006|) . However, despite the diversity of 
the methods, most were applied only for a few binary systems, 
and we do not know that any of the direct methods would have 
been applied for any anomalously slow apsidal motion system. 
So we mainly concentrate here on the most usual, perhaps we 
can say: "classical", although in several cases very uncertain way 
of the detection of a third component in an eclipsing binary sys- 
tem, the method of which is based on the light-time effect (LITE) 
manifesting on the O-C diagram. However, the light-time effect 
arises from the varying distance of the eclipsing system from the 
observer, due to the revolution of the binary around the centre of 
mass of the triple system. Consequently, when this motion takes 
place nearly in the plane of the sky, the distance of the eclipsing 
binary remains almost constant, so the amplitude of the LITE 
(related to a' sin /') might remain under the detection limit. (The 
same can be said about the detection via the systematic radial 
velocity variation.) Unfortunately, this is exactly the situation in 
our ideal case for miscalculating the real apsidal motion period 
from a small fraction of a complete revolution. 

At this point we have to stop for a s hort interplay, and note , 
that in the case of AS Camelopardalis, iKozvreva et alj (Il999h 
really reported that the O-C curve shows small amplitude pe- 
riodic variations after the removal of the apsidal motion terms. 
They fitted light-time solution, and found a third body of mass 
(m3 - 1.1-1.7 M0) orbiting in an eccentric ie' - 0.5) orbit with 
period P' = 805 d. However, we also carried out the analysis of 
the O-C diagram, but only a ve ry poor, and con sequently, very 
questionable fit was found (see lBorkovitsll2003h . Nevertheless, 
we used our calculated orbital elements as input parameters in 
our numerical studies. 

The question naturally arises of what could be an efficient 
method to detect perturbing bodies (if they exist at all) in such 
configurations. There are several method s listed in the above 
mentioned paper of Pribulla & Rucinskil ( 12006 ') which could 
be satisfactory for this purpose. Instead of repeating them, we 
would like to suggest a further dynamical method, which is 
based on the direct detection of the perturbations of the ter- 
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Fig. 12. Resonance phenomena between the orbital revolution and the stellar rotation. When the stellar rotation rate reaches the 
average orbital revolution rate, the amplitude of the precession of the stellar rotational axes suddenly grows. This figure shows the 
dissipative AS4 (left) and AS3b (right) runs where the rotational angular velocities initially were set to close to the averaged orbital 
one. 



tiary. As one can see from our results, the variation of the ec- 
centricity e of the binary may reach even some 10"^ during 
some decades in our sample system. Such a variation could 
have been detected by the present accuracy of radial velocity 
measurements. This is why we suppose that the different val- 
ues of the eccentricity obtained for AS Cam from the 1969- 
1971 and 1981 observations reflect real variation instead of sim- 
ple observational errors. Note when Malonev et al. 1991 ana- 
lyzed the 1969-197 1 obse rvations of .Hilditch (.1972a.h) and 
iPadaha & Srivastaval (1 1975b . they found that the orbital eccen- 
tricity of the binary was e ^ 0.14, while ' Khaliunin& Kozvreval 
(Il983 h deduced from the December 1981 photometry already 
e « 0.16. Furthermore, according to recent O-C studies (which 
were cited earlier) w ^ 220° - 240° were found. If we assume 
that the third body revolves close to the plane of the sky, the dy- 
namical argument of periastron {g) should have a similar value. 
Then, according to our Eq. (|54] |. Ae is positive, which is also 
in good correspondence with the measured tendency. In conclu- 
sion, we stress, that it is very important to obtain new precise 
light curves, as well as radial velocity measurements from this, 
and other anomalously slow apsidal motion systems. Moreover, 
as we predict the motion of the tertiary in the plane of the sky, 
we can expect positive results mainly from the newest interfero- 
metric equipments in the near future. 

Finally, we also obtained new numerical results on the in- 
teraction of the orbital evolution and stellar rotation in close hi- 
erarchical triple stellar systems containing an eccentric binary. 
The most important fact is that resonances might occur due to 
the variation of the stellar rotational rate during the dissipation- 
driven synchronization process. Such resonances were found for 
example in the case when the rotational rate of one of the stars 
reached the average Keplerian angular velocity of the orbital rev- 
olution. It is necessary to investigate these phenomena further in 
detail. 
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Appendix A: Calculation of higher order perturbative terms 

As mentioned in the text, in the most interesting cases the amplitude of the eccentricity variation on the apse-node time scale 
reaches the unity. Consequently, the eccentricity no longer can be considered as constant on the rhs of the perturbation equations. 
Consequently, first the eccentricity e should be calculated as a function of ^ in an iterative manner. 
The corresponding perturbation equation can be written in the following Taylor-expansion: 



= \ [£]o So -1- {[£ + Bi -1- iVi - £Ai] Sq + [£(Bi - EAi)] Si}o -Ae 

eo 



+ i {[(2 - 2Ai)(Bi +Ni- SAi) + B2+N2- &A2] So 

-I- [2(Bi + Ni)(Bi - EAi) + &{B2 - EA2 + 2Bi - 2EAi - AA^Bi + AEA\)\ Si 

+ \[&{B\ + E^A\-2EAiBi)]l.\ ^{Aef\, (A.2) 

^ h^o ] 

where 

E - \. (A.3) 

6 = \. (A.4, 
1 dA 

A, = (A5) 

= (") 
1 dS 

Bi = e-—, (AJ) 
A de 

, 1 d^B 

M = (A.9) 
A de 

N2 = e (A.10) 

So = (A.11) 



S 



\ + E cos 2g 
dSo 
d£' 
d^So 



d£2 ■ 

(In the calculation of the derivatives above, it should also be considered, that, as it can be seen easily 



(A. 12) 
(A.13) 



' (A.14) 



de \-e^\ C2 
dCi e 

de I- 
dCo 

-r^ = 0, (A. 16) 

de 

as far as the present approximation is used.) As the tidal terms depend on the fifth power of the separation, and, consequently, this 
produces a very strong eccentricity dependence, the Ai, A2 derivatives may produce some numerical problems even in the case 
of medium eccentricities. Nevertheless, in the present situation they have the same order of magnitude as E or £, at least at the 
near-perpendicular configurations. Carrying out the iterative integrations, and writing e into the following Fourier- series 

e = eo{l + eo + ^cos2g -1- £4 cos 4^ -1- eg cos 6^ -1- eg cos 8^) -1- .... (A.17) 

the coefficients up to fifth order in £, £, Ai , A2 are as follows: 

ll\,/3 33 3 \4/9 3\,/9 9\,, 



eo = Al £^H- Al A2 £^-H Ai A2 -h Ai\E''&'- 

" " ' \512 512 128 / I 256 256 / 1128 64 ' 
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11 3 3 3 \ 1 

+ — Bi + —Ni + — B2 + — A^2 S {Bi +Ni)E 

1 32 16 64 64 / 32^ ' 



£ cos 2go 



16 



fi cos 4^0 



5 



31 



256 



1 



1 

32^ 

11 



11 1 1 \ 1 

96 48 192 192 / 32^ ' 



12288 1024 
£ cos 6^0 



3072 



1 \ 3 / 1 29 
,3072 ~ 3072^^' ~ 768"^^)^ ^ \~256 ^ 1536"^' 



Azi ££- + ( i - —E^ 

1536 ^/ \768 128 '/ 64 



& cos 8^0 



1 



61440 



1 



3072 



-E&' + 



3072 



£^£^ E^E> + 

768 



1 

160 



£cos 10^0. 



(A.18) 



2 \64 64 64 V \ 32 32 j 8 1536 384 96 192 16 



1 



1 



—i^Bi +B2 + N2)&' - — (Bi + Ni)E& 



64 



32 



^A, + -A2 -{Bi + NO --(B2 + N2) 



4. 



(A. 19) 



+ — Al - — A2] 8" + (-— Al - — A2] ££3 + [1 _ Aai] £2£2 l^Sg + ^ ^ 

^ 384 192 768 ) \ 192 64 ^) \24 48 j 16 16 

'1 _ - 1a2 + -A^]g' + IaiUs- — £4 + l£2£2 ^ l^Sg 

1I6 4 ' 16 16 '/ \8 4 '/ 384 24 16 



+ T2i6Bi + 4Ni +B2 + N2)& + liBi + Ni)E 
16 8 



4, 



(A.20) 



ee = 



1 

192 
1 



48 



4096 



-£= + 



1024 



192 



-E^G 



64 









(-1 






\ 32 


96 V 



24 



2o3 



1024 



-E'6 



256 



£^£2 £"*£ ■ 



32 



— (8B, + 4A?, +82+ N2)&' - — (Bi + Ari)££ 



192 192 48 / I 32 96 96 / I 24 8 



6a, 



(A.21) 



\3072 3072 768 / \256 1536 512 / \768 384 / 64 

+ ( -^6* + —E6' + —E^6^ + —E'&] 6b, (A.22) 
\3072 256 768 64 j^' ^ ^ 



- -61^^' - 30^^^' - - - T^^'^- 



(A.23) 
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For ^ a similar equation can be written as Eq. (|A.2| i and this gives the u(g) relation as follows: 

n'(M - M,*)) = g + y2 sinlg + 74 sin4g + ye sin 6^ + yg sinSg + yio sin lOg + 

where (D*) ' gives the apsidal motion period in the invariable plane in the unit of the inner orbital period, i.e. 



n^' ■! 1 + I -— A2 + - -AiES + 



1 

64' 



1 

16" 
7 

256" 



_9_ 

128" 



1 

2' 
15 

256" 



1 

1024"'" 512"'' ■ 512"''"'" ■ 2"56 
5 . 45 , 103 

— Al At h 

64 256 128 



+ 1 _A, - —A2 + —A] + ^A,A2\e6' + (--Al - — A2 + —A^Ae^S^ - -A,E'& 



+ + —i-Bi + B2)E? + ( — Bi + —bA E&^ + —ByE^B 
4 128 \64 64 / 16 

-Ai +(--A2 + -A? + -A1A2- -A^'lfi^ + (--Ai - -A2 + 2A?|£fi-Ai£2 + — A2£'* 
\ 8 4 ^ 8 4 V \ 2 2 '/ 512 



--A2+A\ + \-^A2 + -A\ + -A1A2 + ^A;j£^ + [--A2 + 2A\ + -AiA2|££ + [--A2 + -A^E' 



1.2 1 



t2\ 172 



1 1 
+ -B26+ -B2E 

4 ^ 2 ^ 



where as before 
n = AVi -£2. 

Furthermore, 

1 1,1, 
^28 16 



+ 7^1^ - (^1 - ^2 + 2Afj£' + ^ (3Ai + 5A2 - 13A2)££2 + Ia^E^S " i (^Bi + 3^2) fi' - ^BiE& 
^ i/s ^ ' 64 ^ ' 16 64 3z 

i (Ai + A2 - 2A2)£ - -L (Al - A2) £^ + ^ (3Ai + I3A2) £6^ + ^ (a^ + A2) E^& + ^AiE^ 

i(A2-2A2)fi+l(A2-2A?)£-iB2 



-Ibi + Ui(Bi+Ni) 



128 
eo + 



16 
4, 



1 2 1 4 

^^=8^^ 16^ 

-— (Al +A2 +3A?)£2- — Ai£fi+ ^— (Ai + A2)fi'* + — (2Ai - 9A2)££^ - — (8A1 + llA2)£^fi^ - — Ai£^ 
64^ 12 u 22 1 1536"- ' 2; 7ggV 1 2; ^^2"^ 1 2; ^^2 ' 



-^(Ai + 3A2 + 9Al + 3AiA2)£2 - ^(Aj + A2 - 3A?)££ - ^Ai^^ + + B2)& + ^-B,E 



1 



1 



-^A2&^ - —A2E6 - -A2E' 



32 



32 



1.1, 

76 = £■ E^ 

24 32 



1 

1152' 
1 



+ — — (Al + 3A2 - lOA? - 4AiA2)£^ + 7^(9Ai + 9A2 - 35A^)£fi2 + _Li.Ai£2£; _ J—{b^ + B2)&^ - ^B^ES 

576 144 



13 



1 

192' 



96" 



13 



■(Al + 7A2)£" + — (Al + 3A2)££^ + (Ai + A2)£:^£ + -Ai£^' 

1152^ ' 64 ' 144^ ' 8 
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1 

78 = — 



64 



effects in the perturbations of hierarchical triple stellar systems 

and, finally, 

ITul = n*Mo - (go + n sin 2go + 74 sin4go + 76 sin 6^0 + yg sin 8^0 + Tio sin lOgo) + 0{e^ , E^). 
As the next step we carry out the inverse transformation. Introducing the variable 
g - TV(u - Uq), 

the Fourier coefficients of the following equation 

g = ^ + G2 sin 2^ + G4 sin 40 + Gf, sin 60 + sin 80 + Gw sin 10^ + ... 
can be calculated as e.g. 

G„ = - r (g-^)sinH^(^)^dg, 
7!" Jo 

where both sin nQ{g) and ^ 



^ can easily be calculated fr 

1 1,1- 
G2 = + - 

^ -7 C 



(A. 30) 
(A.31) 

(A.32) 
(A.33) 
(A.34) 

(A.35) 



from ] 



+ 



Eq. ( |A.24t . The individual coefficients are as follows: 

8 iO 

-Ay&+ —{2Ai -2A2 +5A?+AiA2)£^ - — (7Ai + IIA2 - 28A?)££2 - — Ai£^£+ — (Bi + 3B2)£^ + — 
4 256 112; ^28^^ ^' 64 64^ ' 16 

-i (Ai + A2 - 2A2)£ - 1ai£ + ^(Ai - A2)£^' - -^(JA, + 29A2)ES^ - ^(A, + A2)E^S - ^ 
- ■ .i(A2-2A2)£-i(A2-2A?)£ + iB2 e^, 



1 1,1, 
-E + -E^ + —E^ 

2 8 If^ 

-iA,3.± 



■BiE& 



-A,E' 



1 1 
+ --A,(B 



+ Ni) 



eo + 



(A.36) 



G4 = + 1^4 

8 16 

.1(A,.A2.A?)£^-| 
1 



1 



1 



-l.Ai(Bi+Ni)& 
64 



1 1 95 

^(14Ai + 15A2)££^ - — (Al +4A2)E^&^ - ^^e'S 



LO 

^(Ai + 3A2 - A? + 5AiA2)£2 - ^(Ai + A2 + 3Af)££ - ^A^E^ - -^(B, + B2)& + ^B^E 

^A2&^ - ^A2E& - Ia2E'- 
32 32 8 



(A.37) 



1.1, 
^^-24^-^32^ 

"2^^^^' ^ ^ ^ 19AiA2)fi^ + Y]^(9Ai + 9A2 
-y^(A..7A2)£^.^ 



-4(Ai + 3A2)££2 - ||(Ai + A2)£'£ - 
128 j/6 8 



1 ^ 



(A.38) 



+ 



1^ 

— —{Al + 7A2)£** ^—{Ai + 3A2)E&^ + — 

24576 18432^ ' ^ 614' 



23 , - 97 , 

(A, + A2)£^£^ AiE^&, 

6U4 ' 4608 



Gio - — 



^5 

Eq. ( IA.34b the other orbital elements can also be easily expressed as a function of u, namely 
'00 + E2 cos 20 + E4 cos 40 + E(, cos 60 + E^ cos 80 + Eio cos 10^) + 



(A.39) 
(A.40) 

(A.41) 
(A.42) 
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where 

£00 = -^Ai£2 + IeS + ^(Ai - A2)S' - ^(2 + 14Ai + 13A2)EG' + ^(4 - 37A,)E^6^ + ^E^S 



1 



1 



1 



1 



1 



1 - -(2Ai + A2 - 3A2)£2 + 1(1 _ 2Ai)££ - ^ES' + ^E^S" + ^E'& + -B,& + -{B, + Ni)E 



1 



1 



--(Al + 2A2)£^ - -(2Ai + A2)E£ 
o 4 



4, 
1 



= - 2^ ^ 128 ' - '^^^^ - 64^^ - '^^^^^ - 8^ ^ 1536 
lo 04 32 



3o2 



192 



384 



+ 



1 



1 



1 



--(1 - Ai)£ + —(1 + 2Ai + 4A2)£^ - —(6 - 32Ai - 7A2)££" - -(1 - 3Ai)£^£ - -(fii + NO 



_(2Ai + A2 - 2A^)£ -^(Bi+ Ni) - -(B2 + N2) 



12 



(A.43) 



(A.44) 



E4 = ^(i+a,)S'-Ies- 

lo 



^(6 + I6A1 - 3A2)£^ + ^(3 + 2Ai + 6A2)££^ + ^(3 + 3SA0E^S' - ^E^S 



1(1 . A2 - 3A?)£^ - 1(1 - 2A,)i:£ - 1^£^ . ^E8^ . ^-^eW - j-^E^S 



+ l(fii+M)£ + 

+l(2Bi + 4Nr +B2 + N2) - i(Bi + Nr)E 
lo o 



eo + 



l(-2Ai + 5A2)£' + i(2Ai + A2)££ 



(A.45) 



Ee = -3^(2 + 7Ai + A2 + 3A2)£3 + J^(6 + lAOES" - ^E^S + 
-ls'£ - i^(2Bi + +B2+ N2)e? + l(Bi + Ni)E& 



—£5 - —££4 + ^£2£3 + ii£3g2 

4096 3072 1024 768 



+ 



-— (1 + 4Ai + 8A2)£^ + —(6 - 4Ai + 7A2)££^ - 1(1 - 3Ai)E^& 



(A.46) 



3072 256 2304^ 64 



(A.47) 



" 61440 3072 3072 768 160 



(A.48) 



The perturbations of the other orbital elements can be calculated in a similar way, but as our final purpose is to obtain the 
analytical form of the O-C in the function of the eclipsing cycle number (which is highly related to u, as well as 0), now we use the 
following direct relation, 



dX (dX 



1 ( dX 



dg [dgj^^ ^ (de d^),„ ^ 2 [de^ dg]„, ^^'^ 



(A.49) 



where X means any of the remaining orbital elements or related quantities. So, for the node (h) in the dynamical system, as well as 
for Mmi (which the latter denotes that part of Mm which can be derived from dh cos ii): 

h = h*Q + Ho0 + H2 sin 20 + sin 4^ + sin 60 + //g sin 8^ + (A.50) 

= (Mnii)o + f/o^+f/2sin2^+f/4sin4^+f/6sin6^+f/8sin8^ + ..., (A.51) 

where 



He. 



J + K + 



^{K2 + J2)-\a,{J,+K,) 



£2 + ^{Ji+ Ki)E& 



+ ^(-AiK + KOS + ^KE - -^{K, - K2)e? + ^(Ki + 2K2)ES' + ^K^E^S + ^KE^ 

J1+K1 + h^Ki +K2- 2AiKi - AiK - A2K)& + ^{Ki - AiK)E + ^(^2 + K2)&^ + ^{Ji +K1+J2 + K2)E& 



\{J2 + K2) + ^K2& + ^-K2E 



el 



(A.52) 
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- +K1+J2 + K2- ArJ, - ArKi}6 - ^(K, + 1K2)&' - ^{K, + K2)E& + \k,E^ 
2 4 04 32 6 



-\k2-\{J2 + K2)& 



4, 



(A.53) 



Ha = —{Kx + AiK)& - -KE - —{Ki - 3X2)8^ + —(2Ki + K2)E&^ - —KiE^S 
16^ ^ 8 768^ ^ 128^ ^ 192 

+ +Kx+J2 + K2+AiJi+ AiKOS" - + Kx)ES 

04 32 



+ 



1 



1 



1 



^ ^(z:! + /i:2 + Al/: + A2/:)£ - -(Ki - AiK)e + —iJi+Ki+ 3J2 + 3K2)&^ 

16 8 64 



32 



(Ji +K1+J2+ K2)E& 



^K2e-^K2E 



4> 



(A.54) 



He ^ - 



— (Al + A2)K + —(1 + 3Ai)Ki + —K2 
384^ ^ 192^ ^ 192 ^ 



£2 + — (7Ai/: + 6Ki)E& - —KE^ 
192^ ' 24 



Y^iJi +Ki+ 3J2 + 3K2)6^ + ^(/i + /STi + 72 + K2)E8' - ^(7i + /sT^E^S 



1 



1 



1 



'192^'^' ^ ^^'^ ^ 32^^^' ^ '^^^^^ " 24^'^^' 



(A.55) 



where 



(A.56) 



^2 



„ 1 dAhi 
'n* de ' 
„2 1 d'Ahl 

1 



^ = H^A,2, 



/sTi = e 

and 
Ahi 



1 dAh2 

n* de ' 
2 1 d^Ah2 



COSJi 

2 



-A„i 



while Ah2 could be derived from An2 in a similar maimer. Here J was treated as second order in e, while 7i_2, 
considered as third order. To obtain the corresponding expressions for Wmi, Ahi,2 and its derivatives should be 
-Ani,2 and derivatives. As these quantities will be used later, for the sake of the clarity we define them here: 



(A.57) 
(A.58) 
(A.59) 
(A.60) 
(A.61) 
(A.62) 



L2 



„ 1 dAni 

'n* de ' 

„2 1 d^Anl 

° n* d^e ' 



(A.63) 

K, and Ki^2-s were 
simply replaced by 



(A.64) 
(A.65) 

(A.66) 
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Ml = -e 
Ml = -e 



1 dA, 



n2 



n* de ' 

2 1 d2A„2 
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(A.67) 
(A.68) 
(A.69) 



Formally, similar expression can be written for the direct perturbative terms in the orbital motion (5). Nevertheless, in this case 
the magnitude of the derivatives differ from those above, so we rewrite the results according to the orders of the direct terms as 
follows: 



S = 6q+ Dog + D2 sin 2^ + D4 sin 4^ + De sin 6^ + sin 8^ ■ 
where 



Da 



1 



1 



V + \—V2- -AiVi\&' + -ViE& 



1 



+ 1 + AiWo)£ - ):WE + -i-CWi - W2)£? - ^(Wi + 2W2)E&^ - ^W^E^e, - \wE^ 

4 Z iZo jZ dL o 



1 



1 



Vi + ^(^2 - lAiYx - A2V1 - 2A1V2W + 7 (V + 1 + V2 - 2AiVi) E& 
o 4 



-^{Wi+W2-AiW-A2W-2AiWi)£,- h^Wi-AiW)E 



11,1 1 1 

-V2 + — ViS + -V2E& - -W2& - -W2E 
2 16 4 4 4 



(A.70) 



(A.71) 



D2 



-^Vifi + ^(2^1 - 2y2 + SAiVi + 3A2V1 + 6AiV2)£? - ^^(6^1 + 6V2 - 7AiVi)E6^ - j^ViE^& 

+ + -7^(2^1 + 6W2 - A1W-A2W - 10AiWi)£^ + J-(10Wi + 7AiW)E6 - ^WE^ 
2 128 64 8 

+ -jCVi + V2 + Aivos + - y2)£3 - A(yi + 3y2)££' - ^(Vi + y2)£'£ - iyi(Bi + n^) 

4 iZo 04 lO 4 



Z 04 JZ 



63 



-^(2y2 - 2Aiyi - A2yi - 2Aiy2)£ + ^W2 

o 4 



4> 



(A.72) 



^4 = ^(Vi + V2+ AiVi)3^ - ^ViE8 - -^(V, + V2)G' + :^{Vi + 2V2)E6' + ^(Vi + V2)E^£^ - ^V,E^& 
64 JZ 1j36 Z56 Jo4 64 

+ 7^7 + iVi)e - 4(^1 +MW)£> + \WE + - 3^2)6^ - t^(2Wi + W2)E&' + t^tWi^^S 

o4 lo o /DO IZo lyz 



1 



1 



1 



1 



— (yi + 3y2+A2yi -Aiy2)£''- — (yi + y2-2Aiyi)££- — (Wi + W2 + A1W+A2WOS+ -(Wi -AiW)e 



32 



16^ 



8 



1 



32^26^- ^V2££- ^^26+ ^^2£ 



63 



(A.73) 



^ (2yi + 6y2 + 7Aiyi + lAiyi + 6Aiy2)£^ + y^^^^i + 6^2 + lAiVi)E&^ - Y^ViE^B 



2304 

+ T^(2Wi +2W2 +AiW + A2Vy + 6AiWi)62 - --i-(6V^i + 7AiW)££ + -^E^^ 

"Tlk^^' ^ ^^^^^^ ^ 1^2^^' ^ 3y2)££2 - Y^(yi + y2)£2£ 

+ Y^(Wi + 3W2)£2 - ^(1^1 + W2)££ + Ya^xE" 60, 



(A.74) 



. 7y2)£^ - ^(y. . 3y2)E£3 . ^(y, . y2)E^£^ - ^^v.e^& 



(A.75) 
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In the above equations 

y = ^r^d, (A.76) 

n* 

W = i^Bd, (A.79) 

Wi = e^^^, (A.80) 
n* ae 

where in our computations Vs were considered as first, while Ws as second order quantities in E (or £, e). 

As a next step we calculate the angular elements in the observational system of references. We apply the following relations 
from the theory of spherical geometry: 

cos i = cos lo cos j'l - sin /q sin ii cos h, (A.82) 

„ cos /j -cos /q cos J , sin /() sin /z cos J 

dwm = dh J-. dh — ■ — , (A.83) 

1 - cos^ I 1 - cos^ I 

,. cos/n - cosj'i cos/ , smlnsinh 
dQ. = dh ^ + d/, (A.84) 

1 - cos^ ; 1 - cos^ I 

, cos /n sin ii + cos ii sin In cos h sin h sin /n sin ii 

di = dii : dh . (A.85) 

sm I sm I 

By the use of these relations, after some algebra and the Taylorian expansion of the 1 - cos^ i denominator we obtain 

- dQ cos i = dh cos ji(Coo + C02 cos 2ii + Cqa cos 4/i + C06 cos 6/1) 

+d/icos/isinji(Cio + Ci2Cos2ji +C14 cos4ji +Ci6Cos6ji) 
+d/icos2/zcos;i(C2o + C22COS2/1 + C24Cos4;i + C26COs6ii) 
+d/!COs3/zsin/i(C3o + C32COS2/1 +C34COS4/1 +C36Cos6;i) 
+dh cosAh COS ii(C4o + C42 cos 2ii + C44 cos4ii + C46 cos 6/1) 
+dh cos 5h sin i\ (C50 + C52 cos 2ii + C54 cos Aii + C56 cos 
+dh cos 6/1 cos ;i(C6o + C^i cos 2ii + C(a cos4ii + Cgg cos 

-dii sin/icosji(>Sio +>Si2Cos2ji +5i4Cos4ji) + (A.86) 

where 

765 4725 ^ 357 ^, 91 ^, 

Coo = cos 2/0 H COS 4/0 COS 6/0, (A.87) 

4096 8192 " 4096 ° 8192 ° ^ ' 

1271 8295 1743 329 

C02 1- COS 2/0 COS 4/0 + COS 6/0, (A.88) 

" 8192 16384 8192 " 16384 " ^ ' 

117 525 ^ 483 189 

Co4 = + COS 2/0 + COS 47o cos 6/0, (A.89) 

4096 8192 4096 8192 ^ ' 

25 105 63 231 

= 8192 ^ 16384 ^ 8192 '^"^^'^ ^ 16384 '""'^''^ ^^-^'^ 

3449 31 5 

C12 = sin 2/0 sin 4/o sin 6/0, (A.9 1) 

4096 1024 " 4096 ' 

31 27 3 

Ci4 = -— sin 2/0 - — sin 4/0 - — sin 6/0, (A.92) 
15 9 99 

C16 = sin 2/0 sin 4/0 sin 6/0, (A.93) 

^ 4096 1024 " 4096 

1271 4549 1121 235 

= 8192 - 16384 ^ 8192 ' 16384 ^^"^"^ 

1893 10767 ^ 4003 „ 1025 ^, 

C22 = + cos 2/0 COS 4/0 H COS 6/0, (A.95) 

16384 32768 16384 ° 32768 ° ^ 

287 707 ^ 903 ^ 525 , 

C24 = cos 2/0 H COS 4/0 COS 6/0, (A. 96) 

8192 16384 " 8192 " 16384 " ^ ' 
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75 255 45 495 

C26 cos 2/0 COS 4/0 H COS 6/0, (A. 97) 

16384 32768 16384 32768 

501 57 15 

= -8192 ^^"2''' ^ 2048 ^ 8192 ^^''^'^ 

- ^ sm2/o - sm4/o ^ sm6/o, (A.99) 

45 15 55 

C36 = sm2/„ + — sm4/o - — sm6/o, (A.lOO) 

117 387 99 45 

^ 40% - 8192 4096 '^"^^'"-8192'^°^^'"' ^^'^^'^ 

287 1057 329 175 

C42 — 1 cos 2/(1 cos 4/0 H COS 6/0, (A. 102) 

8192 16384 8192 16384 

19 149 85 59 

^ 4096 - 8192 '''' ^ 4096 '''' ' 8192 '''' '^^^ 

15 15 39 33 

"^'^ = 8192 ^ 16384'^°^ ''"-8192'^°^^'° -^ 16384'^^ ^^''^^^ 

25 5 5 

C52 sin 2/0 H sin 4/o sin 6/0, 

8192 2048 8192 

5 1 1 

C=i4 - sin2/o sin4/o + sin6/(), (A.105) 

1024 " 256 1024 " ^ ' 

15 3 3 

C% sin 2/0 + sin4/n sin 6/0, (A. 106) 

8192 2048 8192 

25 75 15 5 

Ceo = cos 2/() H cos 4/o cos 6/0, (A. 107) 

8192 16384 8192 16384 

75 225 45 15 

Cfi2 + cos 2/0 cos 4/0 + cos 6/0, (A. 108) 

^ 16382 32768 16384 32768 

15 45 9 3 

^ 8192" 16384 8192'^°^"'° - 16384'^°^^ ^^'^^'^^ 

5 15 3 1 

Cfifi = 1 cos 2/0 cos 4/0 H COS 6/0, (A. 110) 

16384 32768 16384 32768 

1727 17 19 

= "2048 ^^"''° -512^^"'''° " 2048 ^^"^'°' ^^■'^'^ 
31 27 3 

S12 — sin2/() sin4/o sin6/o, (A. 112) 

512 128 512 ^ ^ 

5 3 33 

S u — sin 2/0 sin 4/o sin 6/0 (A. 113) 

2048 512 2048 

As Eq. ( |29l ) reveals there is only a very small cyclic variation in /]. Consequently, in Eq. ( IA.86I 1 ii can be considered as constant. 
With this approximation the integration is a very easy task, and we do not feel the necessity to report its result here. (Nevertheless, 
it can be read directly - with the exclusion of the secular term - from the Vonm coefficients [Eqs. (lA. 21 1I )-( IA. 2251 )1 of the final form 
of the O-C given below [Eq. ( lA. 1231 )1.) Instead we list the analytical form of the argument of periastron in the observational frame 
of reference (cS): 

W = g + Um, 

- ojq + Oo0 + O2 sin 20 + O4 sin 40 + Of, sin 6@ + Og sin 8^ 

+0„o sin(n/!(*) + riH^Q) + On+m sin[nh*Q + (nHo + m)@], (A. 1 14) 

where 

(9o = l + f/o(l+Co), (A. 115) 

0„ = G„ + Un{l +Co), (A. 116) 

= -C„, (A. 117) 

n 



e, + lK24\, (A. 118) 



^Ki + +J2+K1+ K2)& 

On^A = ±Cn \~{K, + AiK)S + ^KE - +72+^^1+ ^2)^' + ^(^i + Ki)ES 

I 32 16 1/0 64 
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+ -^iKi+K2)&+^KiE eoj, 
On^6 = ±Cn -^{Ki + K2)&' - ^K^ES + ^KE^ , 



(A.119) 



(A. 120) 



and, furthermore, we applied the following abbreviations: 



Co = Coo + Co2 cos 2ii + Co4 cos 4ii + Coe cos 6ii, 



(A.121) 



but e.g. 



Ci = sinji(Cio + Ci2Cos2ji +Ci4COs4ji +Ci6COs6ji), 



(A. 122) 



i.e. at the trigonometric terms the sin ii or cos ii multiplicators outside the parenthesis are also included. (Note, in the following for 
the sake of the simplicity we consider C„s as being in the order of e". This is not necessarily exactiy correct, but as C„s come from 
the Taylorian expansion of 1 - cos^ i it may be partly verified.) 

We are now in the position to give the final analytic form of the O-C curve up to the fifth order. This is given as follows: 



+y400 sin(4w; + 40o&) + V500 cos(5w* + 5Oo0) 

+Vm cos[w; + (Co + 2)0] + Vi-20 cos[a>; + (Oo - 2)0] 

+yi4ocos[w; + (Oo + 4)g] + Vi-40 cos[w; + (Oq - A)g] 

+Vm) cos[wo + (Oo + 6)^] + Vi-bo cos[W(, + (Oo - 6)^] 

+yi8ocos[wo + (Oo + 8)^] + Vi-so cos[wo + (Oo - 8)^] 

+ V220 sin[2^ + (20o + 2)0] + ^2-20 sin[2wS + (20o - 2)0] 

+V2A0 sin[2wS + (20o + 4)^] + ^2-40 sin[2wS + (20o - A)g] 

+V26Q sin[2w,* + (20o + 6)^] + Va-eo sin[2w,* + (20o - 6)^] 

+y32o sin[3wS + (30o + 2)0] + V2-2Q sin[3w5 + (30o - 2)0] 

+Vim sin[3wS + (30o + A)0] + V2-W sin[3wS + (30o - A)0] 

+V420 sin[4w; + (40o + 2)0] + Vi-2Q sin[4w; + (40o - 2)0] 

+V101 cos[w; +hl + (Oo + Hq)0] + Vio-i cos[w* - /i*, + (Oo - Hq)0] 

+Vw2 cos[w* + 2hl + (Oo + 2Ha)0] + Vm-i cos[w,*, - 2/z; + (Oo - 2Ha)0] 

+Vm cos[w* + 3hl + (Oq + ^Hq)0] + Vio-s cos[w* - 3/zS + (Oq - ^Hq)0] 

+yi04 cos[wo + 4/z* + (Oo + AHo)0] + V10-4 cosLwg - Ah^ + (Oo - AH(,)0] 

+V121 cos[w5 +hl + {Oo + 2 + Ho)0] + Vu-i cos[wS -K + iOo + 2- Ho)0] 

+yi_2i cos[w; + hi + (Oo - 2 + Hq)0] + Vi^2-i cos[w,* -h; + (Oo-2- Ho)0] 

+V122 cos[w5 + 2hQ + (Oo + 2 + 2Ho)0] + Vn-i cos[c4 - 2hl + (Oo + 2 - 2i/o)^] 

+yi_22 cos[wS + 2hl + (Oo - 2 + 2Hq)0] + yi_2-2 cos[w; - 2hl + (Oo - 2 - 2Ho)0] 

+Vm cos[wS + 2>hl + (Oo + 2 + 3Ho)0] + yi2-3 cos[aj; - 3hl + (Oo + 2 - 3Ho)0] 

+yi_23 cos[w; + 2>hl + (Oo - 2 + 3Hq)0] + V,-2-i cos[w; - 3hl + (Oo - 2 - 2,Hq)0] 

+VxA\ cos[w,* + /z; + (Oo + 4 + Hq)0\ + yi4-i cos[w,; - /j^ + (Oo + 4 - Hq)0] 

+V1-41 cos[w; + /i; + (Oo - 4 + Ho)0] + yi_2-i cos[w; - /i^ + (Oo - 4 - Hq)0] 

+Vu2 cos[w* + 2hl + (Oo + 4 + 2Hq)0] + V^-i cos[^ - 2hl + (Oo + 4 - 2Hq)0] 

+Vi-i2 cos[wo + 2/?S + (Co - 4 + 2Hq)0] + Vi-2-2 cos[w; - 2hl + (Co - 4 - 2Hq)0] 

+Viei cos[wo + /iS + (Oo + 6 + Hq)0] + yje-i cos[wS - /z^ + (Oo + 6 - Ho)0] 

+yi-6i cos[w; +hl + (0a-6 + Hq)0] + Vi-e-i cos[w; - + (Oo - 6 - i/o)^] 

+y2oi sin[2w; + h; + (20o + //o)^] + V20-1 sin[2w; -h; + (20o - Ho)0] 

+V202 sin[2w; + 2/j; + (20o + 2Ho)0] + V20-1 sin[2a>; - 2hl + (20o - 2Ho)0] 

+V203 sin[2w; + 3hl + (20o + 3Ho)0] + V20-1 sm[2col - 3h; + (20o - 3Ho)^] 

+y22i sin[2w; + hl + (20o + 2 + Ha)0] + V22-1 sin[2w* -h; + (20o + 2 - Ho)^] 

+y2-2i sin[2w; + /zj + (20o - 2 + //q)^] + V2-2-1 sin[2wS - /z^ + (20o - 2 - Ho)0] 

+V222 sin[2w; + 2/i; + (20o + 2 + 2Hq)0] + y22-2 sin[2a>o - 2hl + (20o + 2 - 2//o)6'] 

+y2-22 sin[2wS + 2hl + (20o - 2 + 2Hq)0] + V2-2-2 sin[2wS - 2hl + (20o - 2 - 2//o)^] 

+y24i sin[2wS + /zS + (20o + 4 + Ho)^] + y24-i sin[2wS - ^^S + (2Co + 4 - //o)^] 

+y2-4i sin[2wS + hl + (20o - 4 + 7/o)^] + V2-4-1 sin[2wS - /zj + (20o - 4 - Hq)0] 
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+y3oi cos[3co; + h; + ooo + Ho)g] + Vjo-i cosp^; -h; + (30o - Ho)g] 

+V302 cos[3w*, + 2hl + (30q + 2Ho)g] + V30-2 cos[3wS - 2h; + (30q - 2//o)^] 

+Vm cos[3w; + h; + (30o + 2 + HqW] + V32-1 cos[3a;; -h*Q + (30o + 2- H(,)g] 

+¥3-21 cos[3w; + hl + (3(9o - 2 + HQ)g] + Vj-z-i cos[3a>* -hl + (3(9o - 2 - //q)^] 

+Vm sin[4w; + hl + (40o + //o)^] + V40-1 sin[4w; -hl + (40o - //o)^] 

+yo2o sin 2g + yo4o sin Ag + yo6o sin 6^ + Vogo sin 8^ 

+ yooi sin(/!*, + //o^) + yoo2 sin(2/!*, + 2//,,^) + Vooa sin(3/i* + 3ffo^) 

+yoo4 sin(4/z; + AH^g) + Vqqs sin(5/z; + SH^g) 

+Vmi sin[/i; + {Hq + 2)g] + Vo-21 sin[/j; + (Hq - 2)g] 

+V022 sin[2/j; + (2//o + 2)g] + Vo-22 sm[2h*o + (2//o - 2)^] + 

+yo4i sm[h'o + (Ho + A)g] + Vo-41 sm[hl + (Hq - A)g] + 0{e\ E\&\ 



jeo 1-2 + — (8A1 + a2)£2 _ 1(4 + Ai)££ + -E^ + —ES' - —E^S" - —E'& + —E^ 
J u| 32*- i u u g g4 256 ^'^'^ 



128 



512 



'1 1 



1 



+ 1 - + -E& - —E^Cj - —CI + -Ci - -ME{\ + Co) 



2 8 



32 



1 ^ . 1^2 1 



32 



8 



-2 + i(2Ai + A2)E^ - 1(4 - 6A1 + A2)E6 + 1(1 - 2Ai)E2 + Ic^ 



jeo ( 1(2 - Ai)£ + Ie - -l-(4 + 7Ai + 9A2)£^ + -l-(20 - 71Ai - 21A2)££' + -^(15 - 22Ai)£2£ 
A I z5o ZDO 04 



^(2-A0£.iE 



1m+1(Li+Mi)£ 



(1+Co) 



1(2 - 4A, - A2 + 2A?)£ + 1(1 - A\)E - —£? + —E£} + —E^E, + —E^ + B, + 1^ 
4^ 2^ 64 64 64 64 2 



1 „ 1 



1 



-£ + -£ Cf - -(M + Mi)(l + Co) 



1 



eo ■ 



-l(Ai+A2)£-l(2Ai+A2)£ 



1 



1 



jeo 1 7(2 + Ai)£ - Ifi - —(4 + 13A| + 3A2)£' + —(28 + 25Ai + 23A2)££' - —(1 - 17Ai)£^£ 
4 2 256 256 64 



-— (2 + Ai)£+ Ifi 
16^ ^' 8 



1m+1(Li +Mi)fi 



(1+Co) 



1 



1 



-(2 + A2-2A0£-;r(l-Ai)£- 



64 



££^ 



64 



64 



1 „ 1 



1 



l(-Ai+A2)£+l(2Ai+A2)£ 



=0 r ' 
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J.o \ ^(-4 . A, . A2)£^ - -Am . ^i^^ . - - ^^E^E.^ . ^^i^^^ . 



J_o2 _ J_c2l/^2 

64 64 



£^|Cf + |lM+^Mi|£(l+Co) 



±(_2 + 5Ai + 2A2)£2 - ^(2Ai + A2)££ + -^(1 - 2Ai)£2 

JZ JZ lo 



1 



1 



19 



m\ (4 + 9Ai + A2 + 2Af )£^ + —(8 + 5Ai)E£ E^ + & E& + 

64^ ^ ^ ^' 32 ^' 16 384 48 768 



£2£2 + 



192 



1 



16 



1 ^2 1 



_£Ar, + (_£^-_££+_£^|Ct + 



16 



64 



1 



1 



-M+—Mx\&+-ME 
8 16 



1 



(1+Co) 



~{2 + 5Ai + 6A2)£2 + ^(8 - 6A1 + 5A2)E£ + ^(-\ + 2A{)E' 
il JZ io 



;eo 1^ (12 + 13Ai - 9A2)£2 + ^(-12 + 5Ai + 3A2)£'£ - ^(3 + 7Ai)££2 + ^E^ 
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Vi4i = jeoCi {- Y^(4 + 9Ai + A2)&^ + ^(8 + 5Ai)££ - + |-^£2 ^ 1 ^ ^£2j ^ ^^ ^42) 

= ^oC. . ^^S^E . ±-^8E - ^^Ej , (A.143) 

Vr-. = jeoC, - ^^8^E - ^^8E . ^^^} , (A.144) 

f 1 3 23 5 1 

V.-. = jeoC, |-^£^ . :^£^E - ^££^ . , (A.145) 

Vio-2 = jeo ||^C2 - \ctj |l + 60 + - ^£ j + ^C.Cs - ^C?C2 + ^Ct} , (A.147) 

Vio2 = jeo |(^C2 + ^C?J (-1 - eo - + ^^'j + ^CiCs + ^C?C2 + ^Clj , (A.148) 

Vi-2-2 = jeoiCl - 2C2) 1^(2 - Ai)£ + + (^S + ^^j foj , (A.149) 

Vi-22 = jeoiCl + 2C2) {^(2 - Ai)£ + ^£ + + eoj , (A. 150) 

V12-2 = jeoiCl - 2C2) 1^(2 + Ai)£ -^E + |^£ - ^isj eoj , (A. 151) 

V122 = jeoiCl + 2C2) 1^(2 + Ai)£ -^E + ^^£ - ^Zij eoj , (A. 152) 

Vi.4-2 = yeoCC? - 2C2) {"Y^S' + Y^^'} ' (A.153) 

Vi-42 = >o(C? + 2C2) {-^S' + Y^£'} , (A.154) 

V14-2 = jeoiCi - 2C2) {-y^e' + - > (A. 155) 

Vi42 = jeoiCi + 2C2) {-y^fi' + - ^£^1 , (A.156) 

yio_3 = ^J^oiC] - 6C1C2 + 8C3)(1 + eo), (A. 157) 

V103 = -^JeoiCl + 6C1C2 + 8C3)(1 + 6b), (A.158) 

Vi-2-3 = —jeoiCl - 6C1C2 + 8C3)(£ + E), (A. 159) 
yo 

Vi-23 = ^jeoiCl + 6CiC2 + SC3)iS + E), (A. 160) 

^12-3 = -ri^eoCC? - 6C1C2 + 8C3)(£ - £), (A.161) 

V123 = ^jeoiCl+6CiC2 + SC3)iS-E), (A. 162) 
yo 

Vio-4 = -^jeoi-Cl + 12C?C2 - 32C1C3 - 12C^ + 48C4), (A. 163) 

Vm = - -T^jeoiCl + 12C?C2 + 32C1C3 + 12C^ + 48C4), (A. 164) 
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V200 = ^e2|l + ie2 + l(l-2Ai)£2 + i(2 + Ai)££-i£2_c2 



2+jel+^6' + E6-^E^-2Cj 



^2-40 
V240 



V: 



V2-20 = el + 3Ai - 2el)S - ^(6 + el)E - ^S^E - ^^6E' - ^E' 



1(_4 + + 3A2)£ + ^(-2 + A,)E 



V220 = ^0 i -^(6 + 3Ai + 2e2)£ + 1(6 + el)E - ^£?E + ^&E^ 



-4(4 + Ai+A2)£+|(2-Ai)£ 
Id o 



^0(2!^ 



3 „2 3 



en < —(8 - 5Ai - A2)£^ + 128 " ^'^^^ ^ (16 ^ 16^^ ' ^ ' ' 
eo ] J^(8 + 13Ai + A2)£' - -^(12 + 1A,)&E + ^E^ + {^^6^ - ^&E +Ie^U\ , 
1 



2-60 - 



" I 64 64 



^20-1 
^201 



V2-2-1 
V2-21 
V22-I 
V2-21 



lelCi 



1 - -Ai|£ + £: + 2(£ + £)^ 



-77^0^1 



1 - -Ai|£ + £ + 2(£ + £)eb 



1 + -Ai|£-£ + 2(£-£)eo 



1 + -Ai|£-£ + 2(£-£)eb 



V2-4-1 
V2-41 
V24-1 
V241 



-^elCiiS'+GE), 



3_ 
32 



i^elCii&^+&E), 



-—elCx(&-Z&E^lE\ 
^elCi(S'-3SE + 2E\ 



V20-2 = -eliCl-Ci)il+2eo), 

V202 = lel(Cl+Ci)(l+2eo), 



T. Borkovits et al.: Tidal and rotational effects in the perturbations of hierarchical triple stellar systems 31 

V2-2-2 = -^el(C\-C2){& + E), (A.184) 

V2-22 ^ -^el{C\+C2){& + E), (A.185) 

V22-2 = ~el{C\-C2){£'-E\ (A.186) 
lo 

V2-22 = -^el{C\+C2){&-E), (A.187) 

V20-3 = -lel(C\-3CiC2 + 2Ci), (A.188) 

o 

V203 = ^4C^ + 3C1C2 + 2C3), (A.189) 

V300 = \jel |l + \el + ^6' + ^E6 - ^E^ - \c\ + 3eb + S^j , (A.190) 

V3-20 = ;e^|^(-2 + Ai)£-i£ + (-^£-^£jeo|, (A.191) 

V320 = 7e^|-i(2+Ai)£+i£ + |-^£+^£jeb|, (A.192) 

V3-40 = el |^£2 + l££ + Ie^I ^ (^ ^93) 

V340 = - + , (A. 194) 

V30-1 = -iye^CiCl+Seo), (A.195) 

V301 = ^Je^CiCl + 3eo), (A.196) 

(A. 197) 

y3-2-i = ^7£^Ci(£ + £), (A.198) 

y3_2i = -^;£^Ci(£ + £), (A. 199) 

V32-1 = \iElCx{e.-E), (A.200) 

V321 = -^7£'oCi(£-£), (A.201) 

y3o-2 = ^ je^(3C? - 2C2), (A.202) 

V302 = \iel{l>C\ + 2C2), (A.203) 

V400 = -^4(l+4eb), (A.204) 

(A.205) 

V4-20 = ^^{^e+^s}, (A.206) 

V420 = (A.207) 
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4 

5 4 

V401 = "Ye^"^^' 



(A.208) 

(A.209) 



V500 = -^jeo- 



(A.210) 



V020 = -7 Vi£ + ;^^(2Vi - 2V2 + 5AiVi + 3A2V1 + 6AiV2)6^ - T^CeVi + 6V2 - lAiVi)ES' - ^ViE^S^ 

+ l;W+ -^(2Wi +6W2-AiW- A2W - lOAiWi)8^ + ^(lOWi + 7AiW)E6 - IwE^ 

2 128 64 8 

-t(^i + ^2 + Aiyi)£ + -^(Vi - y2)£3 - ^(Vi + 3y2)££2 - rj^CVi + V2)E^& - ^ViiBi + Ni) 



128 



1 1 
--{1^2 - 2AiVi - A2V1 - 2AiV2)6 + -W2 
8 4 



4 



+Co ^ + i(Li + Mi)6 + ^(Mi + 3M2)£' + |-Mi££ - ^ME^ 

, Z 4 04 JZ 6 



^Mi + i(Li + L2 + Ml + M2) 



60 + ^M26^| , 



(A.211) 



Vo40 = ^(Vi + y2 + Aiyi)£2 - ^ViES - T^(Vi + V2)£' + + 2y2)££' + ^(^1 + ^2)£'£' - ^V'i£3£ 

04 3z ioio z5o 3o4 54 

+ ^ + A^i)£ - i^CWi + Al W)£ + \v^E + -i-(Wi - 3W2)£^ - -^(2Wi + W2)££' + T^WiE^g 
o4 lo 0/00 1Z5 lyz 

+ ^(Vi +3y2+A2yi -Aiy2)£2- :^(yi + y2-2Aiyi)E£- -^(Wi + W2+AiW + A2W)£+ ^(Wi -AiW)E 
64 32 16 8 

+ ^ V2£2 - ^ V2££ - ^W2£ + ^lE el 
+Co j-^j^W + AiM)£ + 1m£ - ^(1,1 + L2 + Ml + 

V060 = -5^(2^^! + 6y2 + 7Aiyi + lA2yi + 6Aiy2)£3 + -^{(>Vi + 6y2 + 7Aiyi)££2 - ViE^g 
+;-^(2iyi + 2W2 + Al W + A2W + 6A1 Wi)£^ - -^(6^1 + 7Ai W0££ + 

JOT" J. 

+ -Yi^(^i + 7V2)£' + Y^(yi + 3y2)££2 - :i^(yi + y2)£2£ 



h M2)£^ + 



■2 ■ 1(1,1 +Mi)££ 



-i(Mi + M2)£ + 
16 o 



(A.212) 



24 

Y^(Vi + V2)£^£ 

^(Wi + 3^2)6^ - ^{Wi + W2)E& + ^WiE^ 



192 - ^1*^1 + t^2jfio+ ^wi- 

+Co I Y^(Mi + M2)£2 - ^Mi££ + ^ME^j , 



Vooi = -C 



-Ci, 



Vo-21 =-C,[i/^+i 



i/:+i(//i + /s:i)£+i/i:ieo 



(A.213) 



(A.214) 

(A.215) 
(A.216) 
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Vo2l 


= Ci 


Ik + i(Hi + Ki)6 + 

4 






(A.217) 


Vo-41 


= -C 








(A.218) 


Vo41 




^ 32 ' 16 / 
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